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1. INTRODUCTION 


The computer program discussed in this document, CRP, is intended for 
the graphical presentation of a general dependent variable X that is a 
function of two independent variables, U and V. The required input to the 
program is the variation of the dependent variable with one of the indepen- 
dent variables for various fixed values of the other. The computer program 
is named CRP, and the output is provided by the SD 4060 plotter. Program 
CRP is maintained in the GSFC program library and is suitable for operation 
on the IBM 360 series machines. Versions suitable for operation on UNIVAC 
1108 and IBM 7094 machines may be obtained from the National Space Science 
Data Center. 

The dependent variable X = f(U,V) is displayed as the two orthogonal 
families, X = f 2 QJ) for various V = and X = f 2 (V) for various U = Uj . The 
plot is a type of nomograph and is referred to as a carpet plot. An example 
of the output is shown in Figure 1 where the surface magnetic field declination 
given by the IGRF 1965 field model is shown as a function of geodetic lati- 
tude and longitude for the Northern Hemisphere. The display of a variable, 

X = f(U,V), is a classical problem. Among the usual solutions are tabular 
or graphical displays of X = f OJ) for various values of V and contour plots 
of U = f (V) for various values of X. These techniques suffer either from 
being too cumbersome to use and difficult to interpolate or else from 
requiring many pages of display. The carpet plot technique described in 
this document offers a method whereby data may be displayed in a concise 
fashion that is directly amenable to fast and accurate interpolation with- 
out the requirement of a calculator or slide rule. The accuracy with 
which interpolation may be performed is very much dependent upon what data 
are to be displayed and whether the dependent variable is plotted in a linear 
or logarithmic fashion. In many cases an interpolation accuracy of 2-3 per- 
cent may be obtained, and with care and experience even more accurate inter- 
polation may be performed in some cases. 

Program CRP is an extremely flexible program that offers the user a 
wide variety of options. The dependent variable may be presented in either 
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a linear or a logarithmic manner. Automatic centering of the plot is pro- 
vided in the ordinate direction, and the abscissa is scaled automatically 
for a logarithmic plot. Because of the wide variety of plots that may be 
generated by CRP, flexible labeling and addressing of the orthogonal families 
are included. The grid and the plot may be output on successive frames 
of the SO 4060 output in order that color carpet plots may be generated. 

This capability is most useful when complex plots with fine grids are gen- 
erated. Program CRP has been extensively used at the National Space Science 
Data Center, and samples of the applications have been presented by Teague 
and Vette (1972 and 1974). 

In Section 2 of this document, a description of the carpet plot tech- 
nique is given; and in Section 3, the coordinate system used in the program 
is presented. Section 4 discusses various aspects of the program logic 
and gives detailed documentation of the data card format and NAMELIST input 
parameters. A reference summary of these parameters is presented in Appendix 
A. Section 5 is intended to provide the user with assistance in generating 
usable carpet plots. Guidelines for making initial estimates of the impor- 
tant input parameters are given, and methods for improving existing plots 
are discussed. The JCL required with program CRP is discussed in Section 
6, and the program output is discussed in Section 7. 
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2. THE CARPET PLOT DISPLAY 


As noted in the Introduction, the carpet plot display is a form of 
nomograph (Adams, 1964). An arbitrary dependent variable X is presented 
as a function of two independent variables U and V. In theory, an arbitrary 
number of independent variables may be displayed as X = f (U 1 ,U 2 ,U 3 . . .U^) . 

A practical limit of k = 4 exists, and even for k>2 the plots become physi- 
cally very large. The dependent variable is plotted as a continuous function 
of an independent variable at a variety of specified values of the other, i.e., 

X = f(U,V^) = f j (U) for various V = 

The orthogonal set is plotted in the same manner, i.e., 

X = f (Uj , V) = f 2 (V) for various U = Uj 

Throughout this document, the word 'family' is used to refer to the 
set of curves that is X = fj (U) for various V = or the set that is 
X = f 2 (V) for various U = Uj . A single curve of a family is referred to 
as a 'member.' Note that the word 'orthogonal' is used with the strict 
mathematical meaning of independence and does not mean that the two families 
are plotted perpendicular to each other. Conventionally, the dependent 
variable X is plotted as abscissa, and the independent variables as ordinates. 
It is evident from Figure 1 that there is no fixed ordinate axis. The 
origins of the ordinate axis are sliding although the ordinate scales 
remain fixed. These scales are indicated by bars on the plot shown in 
Figure 1. In Figure 2 , an element of a general carpet plot matrix is shown. 
The two members of one orthogonal family X = f 1 (U) for and V^ +1 and the 
two members of the other family X = £, (V) for Uj and Uj +1 define the four 
matrix points (i,j), (i+l,j), (i,j+l), and (i+l,j+l). In this report, the 
convention is adopted that for an arbitrary matrix point (i,j) the first 
index refers to the independent variable V being held constant at and 
the second index to U remaining constant at Uj . Successive points in the 
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matrix are related by the ordinate scales AU and A V and by the differences 
between the values of the independent variables at the successive points. 

In Figure 2, moving the direction of the arrow from point (i,j) to (i,j+l), 
the vertical coordinate of the latter point Y^j+i is 


f iJ+l 


Y i, 


IVl 


-U. 

J 


AU 


and the coordinate of the point Ci + l>j + l) is 


Y. 


i+l,j+l 


= Y. . 
1>1 


|U. -U.l 
1 J + 1 J 1 

AU 




AV 


In this fashion, the entire matrix is built up. The horizontal coordinate 
of each point is determined by evaluation of 

X. . = f (U. ,V.) 
il l i 

Carpet plots are directly amenable to interpolation. In Figure 2, the 
value of the dependent variable is required at some arbitrary point 
P E 0,y) . The intermediate matrix points (i,y), (x,j), (i+l,y), and (x,j+l) 
may be located with respect to the nearest matrix points by equations of 
the form iy y i 

v _ Y + 1 x " 11 

x,j + l Y i,j+1 Av 

In practice, the locations of these points are determined by simple obser- 
vation using the horizontal grid lines (Figure 1). Other intermediate 
matrix points may be determined on other members such as Uj + ^ and ^ if 
the bounding members for point P are very curved and the members U = U^ and 
V = can be estimated to obtain the intercept P. In practice this pro- 
cess can be performed quickly, without the necessity for calculation, and 
the dependent variable X can be found to an accuracy better than that 
given by linear interpolation. As noted in the Introduction, - an accuracy 
of 2-3 percent may be expected for this interpolation, dependent upon 
the data being plotted. 
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3. THE CARPET PLOT COORDINATE SYSTEM 


Program CRP displays on an SD 4060 plotter a parameter X that is a 
function of two independent variables, U and V, in the form of the two 
orthogonal families 

X = f j CU) for various V = V.^ 
and X = f 2 (V) for various U ~ Uj 

For any given point (U,V) X must be single valued. Input to the program 
may be either one or both orthogonal families. In the former case, the 
other family is determined automatically, and this is the normal mode of 
operation of the program. The basic coordinate system used is the SD 4060 
subjective coordinate system, X , Y . (The reader is referred to an SD 4060 
manual for a complete description of the subjective coordinate system.) 

The maximum vertical extent of the plot defining the horizontal bounding 
grid lines is Y = Y as shown in Figure 3. The value of Y is chosen 
arbitrarily to be 100. The horizontal subjective units for X g are taken as 
those of the input dependent variable X. The origin of subject space in 
absolute SD 4060 rasters is the point (500,200), and the bounding grid lines 
in Figure 3 are defined by this point and the point (3600,2500). In the 
raster coordinate system, the maximum extent of the SD 4060 oscilloscope 
screen is defined by the point (4095,3071). 


The vertical subjective coordinate Y^ is related to the real engineering 

units of U and V by the number of horizontal grid lines N and the vertical 

scale parameters AU and AV. That is, in subject space, the distance between 

successive horizontal grid lines is Y^^/N-l and corresponds in real space 

to changes in the independent variables of AU and AV. In Figure 3, the 

curves X = f 1 (U) for V = V x and X = ^ (V) for U = U r define the plot matrix 

point (1,1) having the minimum value of Y g equal to Y . Then the subject 

coordinate Y for X = f . (Y ) for V = V, are 
s 1 s 1 

(U - U.)Y 
^ i max 


Y = Y + 
s o 


AU 


N-l 
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and the coordinates for the ith line V = 


Y 

s 


(U - UJY 
^ l max 

Y o + AU N-l + 


are 


C v. - V )Y 
^ l i' ] 


AV 


max 

N-l 


( 1 ) 


Similarly for the other orthogonal family, the jth line for U = U. is 
given by the expression 


Y 

s 


= Y + 
o 


cv - vjy 


AV 


max 

N-l 


(U. - UJY 

J 1 1 


max 


AU 


( 2 ) 


For ease of interpolation, it is very desirable that the families of plot 
curves intersect at horizontal grid lines, such as the point (i,j) in 
Figure 3. If the point (1,1) is positioned p horizontal grid line incre- 
ments from the line Y =0, then 

s 

pY 

Y _ max 

o “ N-l 


A general parameter ^ is introduced where 


* 


V. U 


l 


1 

ij AV AU 


and equation (1) becomes 


Y _ max 
s N-l 


U V. 
p + + — i- 

' AU AV 


*11 


= f(U) for V = V. 


and equation (2) becomes 


Y _ max 
s " N-l 


U. V 

ip + _ i. + — 
|F AU AV 




f(V) for U = 


U. 

1 


(3) 


(4) 


( 5 ) 


The condition that the arbitrary matrix point (i,j) lies on a horizontal 
grid line is, from equation (4) or (5), that 


P + *ij “ *n 


I. . 


C6) 
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where I is an arbitrary integer and p has already been given as an integer 
quantity. Specifically for the matrix point (n,ra) having the maximum Y g value 


Y = Y 
s max 

N-l 


the condition is that 


o ♦ *„■ - 


p + dj - lb , - I U ) 

V y nm 1 nm 

The parameter p is input to the program and the maximum and minimum matrix 

points (n,m) and CM) are positioned respectively distance Y q below the 

line Y = Y and above the line Y =0. Then, the number of horizontal 
s max s 

grid lines drawn for the plot is 

N = p + I + 1 ( 8 ) 

r nm 

and the condition that any arbitrary point (i,j) lies on one of these grid 
lines is given by equation C6) . If the matrix point (n,m) does not lie on a 
horizontal grid line, the left-hand side of equation (6) is truncated to ob- 

tain I nm’ i,e -’ 

I m ■ IFIXCP ♦ - +„) ' 9 > 

where IPIX is the FORTRAN Fix function. 


7 





4. PROGRAM INPUT AND LOGIC 

Program CRP consists of 12 subroutines, one BLOCK DATA subprogram and 
one FUNCTION subprogram. These are: MAIN, BLOCK DATA, READ, SUPER, INITIAL, 

GRIDPT, SCALE, LINT, REVERS, FINDF , INTERP, CRVFIT, SPLN2 , and SPLIN (Y) the 
FUNCTION subprogram. A detailed description of the functions performed by 
these routines is not given, but reference will be made to them as the 
program input and logic are discussed. A list of the program CRP is 
given in Appendix B for the IBM 360 series machines. 

There are many input variables for Program CRP, and their functions 
may be categorized as follows: 

1 Carpet plot control parameters 

2 Interpolation parameters 

3 Grid parameters 

4 Labeling and Titling parameters. 

As far as possible these parameters are discussed separately, although in 
some cases a parameter may fit into two categories. Most of the input 
parameters are entered via the NAMELIST HEAD. Default values are given for 
each parameter, and in many cases a given parameter need not be input by 
the user. In addition, for cases in which several similar carpet plots are 
generated, the input parameters can be specified for the first plot only. 

The data deck setup is shown in Figure 4. The initial card specifies 
the number of plots to be generated. The variable name is NP LOT , and the 
input format is 13. NPLOT groups of data follow this card. Each group 
contains a NAMELIST HEAD and a set of member data cards containing the 
data to be plotted. 

4.1 The Input Data Cards 

Under control of the NAMELIST variable RTYPE (Section 4.2.1), sub- 
routine READ reads in one or both of the orthogonal families of curves. 
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If only one family is input, the other is determined by interpolation 
(Section 4.2.2). Program CRP assumes that the first family read in is 
the family X = f(U) for various V = V^. If both families are entered, 
the curves X = f(V) for various U = Uj follow. Each set (or pair of 
sets if both families are read in) of member data cards is preceded by 
a NAMELIST HEAD as indicated by Figure 4. 

The deck setup for the member data cards is shown in Figure 5 and 
a detailed summary of the card formats is given in Table 1. From Figure 
5 it is clear that each member of each orthogonal family is preceded by 
a header card that contains the number of data cards for the member and 
the value of the constant independent variable. Each orthogonal family 
is terminated by a blank card. 

The format of the header card information is shown in Table 1. The 
first 67 columns are not used by Program CRP and are available for com- 
ments. Columns 68-69 contain the right justified integer that specifies 
the number of data cards, N s , following for the member. Columns 75-80 
contain, in F6.m format, the value of the independent variable that is 
constant for the members. 

The format for the data cards is shown in Table 1. One of the N data 
cards is shown. Four pairs of independent -dependent variables are input 
on each card with format 4 (F8. 5,1X,E10. 3) . The last data card may contain 
fewer than four pairs, and if this is the case, the fields may be left blanl 
and will be automatically rejected by Program CRP. 

There are two restrictions on the data cards resulting from dimension 
statements. First, a maximum of 50 pairs of points per member is 
allowed with the result that N s <12. Second, a maximum of 20 and a min- 
imum of two members per orthogonal family are allowed. Each member of an 
orthogonal family must be input in ascending order of the constant inde- 
pendent variable, i.e., in Figure 5, for the first orthogonal family, V = 
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Vj, V 2 . . . , must be in ascending order. However, the order in which 
these members are plotted may be specified by the user (Section 4.2.1). 

The varying independent variable on the data cards must be monotonic. 

4.2 The NAMELIST Input 

Forty-six input variables are entered through the NAMELIST HEAD. With 
the exception of the variable IPR, which is discussed in Section 7, these 
variables are discussed in this section under the four headings listed at 
the beginning of Section 4 viz, (1) Carpet plot control parameters, (2) 
Interpolation parameters, (3) Grid parameters, and (4) Labeling and 
Titling parameters. A summary of the NAMELIST variables and their de- 
fault values is given in Appendix A. 

4.2.1 Carpet Plot Control Parameters 

Eight variables are in the category of carpet plot control 
parameters. These are RTYPE, PO, DELTAU, DELTAV, UPLOT, NUPLOT, VPLOT 
and NVPLOT. 


RTYPE is an integer variable that determines whether one or 
both orthogonal families are input to Program CRP. For RTYPE = 0, the 
family X = f (U) for various is input. For RTYPE = 1, this family 
is followed by the family X= f 2 (V) for various Uj as shown in Figure 5. 
For the latter case, the parameters UMAX, UMIN, VMAX, and VMIN (see 
Section 4.2.3) must be specified. These are not necessary for RTYPE = 0. 
The default value for RTYPE is 0. 

The variable PO is the parameter p discussed in Section 3. 

PO controls the vertical positioning of the carpet plot. PO is equal 
to the number of increments of horizontal grid-line spacing that the 
minimum carpet plot matrix point is positioned above the line Y s = 0, 
and the maximum matrix point is positioned below the line Y s = Y max . 

For instance, in Figure 1 PO has a value of six. PO should be entered as 
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a whole number. The default value is four. As discussed in Section 3, 
the carpet plot matrix is generated from the minimum matrix point that 
is fixed on the p+1 th horizontal grid line. The maximum matrix point 
may not lie on a horizontal grid line, and equation (9) is used to deter- 
mine the variable I . As a result, the maximum matrix point will be 

run 

positioned below Y = Y a distance PO minus a fraction of a hori- 
r s max 

zontal grid- line spacing. 

The two variables DELTAU and DELTAV are the parameters AU 
and AV discussed in Sections 2 and 3. These are the most important 
parameters in the carpet plot generation and some discussion of their 
use is given in Section 5. DELTAU and DELTAV not only determine the 
vertical scale for the plot but also influence the number of horizontal 
grid lines that are drawn. As noted, the first orthogonal family input 
is the family X = f j (U) for various The scale parameter for this 

family, DELTAU, must be positive. A restriction on the input is that 
the V^'s must be monotonically increasing (Section 4.1). If DELTAV is 
a positive quantity, the curves for constant V^ will be plotted in the 
order in which they are input with the lowest value of V^ appearing at 
the bottom of the plot. (See Figure 1 for an example.) If DELTAV is a 
negative quantity, the family is plotted in the reverse order, with the 
highest value of V^ appearing at the bottom of the plot. This is a 
useful facility in generating carpet plots and is discussed further in 
Section 5. If both families of curves are input, DELTAU and DELTAV 
should both be positive. No default values are included for DELTAU 
and DELTAV. 


The arrays UPLOT and VPLOT define the curves that are 
plotted on the carpet plot. NUPLOT and NVPLOT are the lengths of the 
arrays UPLOT and VPLOT, respectively. The maximum length is 20. The 
array UPLOT contains the values of Uj for which the family X = f 2 (V) is 
plotted. Similarly VPLOT defines the values of V^ for which the family 
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x = f (U) is plotted. If the single family X = f x (U) for various 
is input, the array VPLOT may be composed of any or all of the values 
of for which X = f x (U) was input. That is, curves may be input that 
are not to be plotted. This is useful when the interpolated family 
X = (V) for various values of Vj has sharp changes of gradient. The 

input curves that are not plotted will provide grid points for the 
interpolation (Section 4.2.2). If both families of curves are input, 
UPLOT and VPLOT must be the same, respectively, as the U.. and V.^ values 
for which X = ^ (V) and X = f x (U) were input. 

There is a restriction on the UPLOT and VPLOT arrays that 
relates to the input data. Whether the family X = f 2 (V) for various 
U. is determined by interpolation or is read in, each line of constant 
UPLOT (I) must be continuous. For instance in the case of interpolation,, 
if ten lines of X = f l (U) are input, and lines 1 through 5 and 7 through 
10 satisfy UPLOT - UjeU, and line 6 does not, program termination occurs. 
However, if lines 3 through 8 satisfy UPLOT = UjeU, and lines 1, 2, 9, 
and 10 do not, the line UPLOT = Uj is continuous, and a correct plot will 
be obtained. If both families are input. Program CRP is not able to 
plot breaks in the input curves as it simply connects adjacent \points 
in the input arrays. 

Note also that if the orthogonal family X = f 2 (V) for var- 
ious Uj is obtained by interpolation, for each value of UPLOT there must 
be at least two grid points given by the input family X= f x (U) for var- 
ious V^ (see Section 4.2.2). 

4.2.2 Interpolation Parameters 

A single parameter, CHOICE, is in this category. T^is par- 
ameter is applicable to the case RTYPE = 0 only (Section 4.2.1), and 
determines which of two interpolation schemes will be used to determine 
the orthogonal family X = f 2 (V) for the array Uj = UPLOT. CHOICE may be 
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entered as either 'CRVFIT' or 'SPLINE'. The default value for CHOICE 
is the former. CHOICE = 'CRVFIT' selects the interpolation routines 
presented by Akima (1969), and CHOICE = 'SPLINE' selects the routines 
given by Thompson (1970). Both interpolation schemes have the charac- 
teristic that, given a set of grid points, a curve that is aesthetically 
pleasing is forced through these points. Note that the schemes are not 
similar to least squares fits as the resultant curve is constrained to 
pass through the grid points. Two different schemes have been included 
in CRP as the authors have encountered problems with spurious oscil- 
lations in separate applications. However, no cases have been encountered 
in which both schemes have been unsuccessful. No prior statement can be 
made concerning the . suitability of either scheme, and this must be de- 
termined by observation. 

For application on smaller computers for which limited core 
capabilities exist, either a dummy subroutine may be substituted for 
CRVFIT for the case CHOICE = 'SPLINE'; or for the case CHOICE = 'CRVFIT', 
a dummy subroutine and function may be substituted for SPLN2 and SPLN(Y), 
respectively. 

For interpolated lines for which only two matrix points exist, 
the above schemes are not used. In this case either linear or logarith- 
mic interpolation between the matrix points is performed, depending upon 
the value of ABSCA. In addition, when only three matrix points exist, 
and CHOICE = 'SPLINE', a parabolic fit is used to determine the inter- 
polated member without entering the spline-fit routines. 

Substituting the input array V^ associated with input family 
X = f L (U) into equation (5), Program CRP determines an array of the subject 
coordinates for the matrix points of the orthogonal family X = f 2 (V) 
for each value contained in the array UP LOT. These arrays of matrix 
points are passed to the curve-fitting interpolation routines to estimate 
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the complete orthogonal family. Each member of this orthogonal family 
is restricted to an array of maximum length 50. The curve-fitting 
routines estimate Np points between the matrix points where 


N p = IFIX 


50 - N v 
N v - 1 


( 10 ) 


where N v is the length of the input array and N v >_ NVPLOT. 


4.2.3 Grid Parameters 

Ten parameters are in the category of grid parameters: UMAX, 

UMIN, VMAX, VMIN, KOLOR, ABSCA, XO, XF, DELTX, and SUBDIV. 

The first four of these variables are required only for the 
case RTYPE = 1, i.e., when both orthogonal families are input. UMAX 
and UMIN are, respectively, the maximum and minimum values of U input 
with either orthogonal family. VMAX and VMIN are similarly defined 
for the independent variable V. 

The integer parameter KOLOR allows the user to split the 
generated carpet plot into two sucessive frames of SD 4060 output. For 
KOLOR = 1, the carpet plot is displayed on a single frame and for 
KOLOR = 2, the plot is split into two frames. The default value for 
KOLOR is 1. 


The- purpose of this facility is to allow the user to gener- 
ate a finished carpet plot in which the grid and the plot are printed 
in different fashions. An example of a finished product is shown in 
Figure 6 in which the grid has been printed in red and the plot and 
labeling in black. The authors have also generated carpet plots in 
which the grid has been printed using a dot matrix film, and the carpet 
matrix has been printed using regular film. This facility is most use- 
ful when complex plots are generated with many horizontal and vertical 
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grid lines, and confusion between grid lines and the members of the 
orthogonal families may result. For KOLOR = 2, Program CRP generates 
a frame of SD 4060 output containing the grid only, followed by a second 
frame containing the plot, labeling, titling, and scale markers. Since 
these two frames must be registered in the course of the enlarging and 
printing processes, corner registration marks are added to the two 
frames for KOLOR = 2. If this facility is to be used, the user is cau- 
tioned that obtaining correct registration of the two frames requires 
considerable care in both the enlargement and printing processes. 

The variable ABSCA determines the abscissa scale for the 
dependent variable X. If ABSCA = 0, this scale is linear; if ABSCA = 1, 
it is logarithmic (base 10). For the former case, the NAMELIST variables 
X0, XF, DELTX, SUBDIV, and XFMT must be specified. These variables 
are not required for ABSCA = 1. The default value for ABSCA is 1. 

The variables XO, XF, DELTX, and SUBDIV (see Section 4.2.4 
for XFMT) are required only for plots with linear abscissa scales 
(ABSCA = 0) . These parameters are determined automatically for loga- 
rithmic plots. XO and XF are, respectively the minimum and maximum 
abscissa values required for the plot. Their units are those of the 
dependent variable X. No default values are given for these parameters. 
DELTX is the abscissa increment (in units of the dependent variable) at 
which the abscissa is to be labeled. Vertical grid lines and tick 
marks are also drawn at these increments. No default is given for 
DELTX. SUBDIV is the number of subdivisions of the increment DELTX; i.e, 
SUBDIV-1 vertical grid lines will be drawn in any incremental abscissa 
distance DELTX. The default value for SUBDIV is one, hence no vertical 
grid lines are drawn between labeled points. Usually the parameter 
DELTX will be an integer multiple of the abscissa distance XF-XO. 

4.2.4 Labeling and Titling Parameters 

There are 26 parameters in the category of Labeling and 
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Titling parameters. They are TITLEI, Nl } TITLE2, N2, TITLE3, N3, XTITLE, 
NX, LABEL1, NLAB1, XP0S1, YP0S1, THETA1 , LABEL2, NLAB2 , XPOS2, YPOS2, 
THETA2 , ULABEL, VLABEL, XTICK, YTICK, DCHAR, NCHAR, CHFMT , and XFMT. 

The first six of these parameters provide a titling capa- 
bility for the top of the carpet plot. The parameters TITLEi may 
contain arbitrary alphanumeric information of maximum length 60 char- 
acters. The parameters Ni contain the number of characters in the corre- 
sponding TITLEi array. TITLEI, 2, and 3 are, respectively, the first, 
second, and third lines of the title at the head of the plot. Each 
line of the title is automatically centered with respect to the grid. 

The default values for the TITLEi parameters are blanks and 1 for the 
N parameters. 

The two parameters XTITLE and NX define the abscissa title 
drawn below the numeric labeling of the abscissa scale. XTITLE contains 
the alphanumeric title with a maximum length of 60 characters, and NX 
is the number of characters in XTITLE. XTITLE is automatically cen- 
tered with respect to the grid. The defaults are blanks for XTITLE 
and 1 for NX. 


The parameters LABELi, NLABi, XPOSi, YPOSi, and THETAi re- 
late to the alphanumeric labeling of the two orthogonal families of 
curves. Since the position of these labels is arbitrary, either set 
of parameters can refer to either family. The alphanumeric title 
information is entered into LABELi and has a maximum length of 36 
characters. NLABi are the number of characters in the LABELi arrays. 

The defaults for LABELi and NLABi are, respectively, blanks 
and 1. The positions of the first characters of LABELi are defined by 
the parameters XPOSi and YPOSi. XPOSi are the abscissa coordinates 
measured from the left-hand grid margin and are entered in inches. The 
range of XPOSi is -1.25 <_ XPOSi/ (inches) £ 8.0. The defaults for XPOSi 
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are zeroes. YPOSi are the ordinate coordinates measured from the bottom 
grid margin and are entered in inches. The range of YPOSi is -0.4 _< 
YPOSi/ (inches) <5.9. The defaults for YPOSi are zeroes. Values for 
XPOSi and YPOSi are given assuming that the SD 4060 hardcopy output 
measures 11 inches by 14 inches (grid size 8 inches by 5.9 inches). 
Occasionally, because of differing enlargements in the xerox process, 
the plot size may differ from the 11 inches by 14 inches nominal, and 
the values of XPOSi and YPOSi must be scaled accordingly. The var- 
iables THETAi define the orientations of the alphanumeric information 
contained in LABELi. THETAi are entered in degrees and define the 
angle measured from the horizontal in a counterclockwise direction. 

The ranges for THETAi are 0 < THETAi/ (degrees) £ 90, and the default 
values are zeroes. 

The variables ULABEL and VLABEL relate to the numeric 
labeling of the members of each orthogonal family. These members may 
be labeled at one or both ends of each curve. VLABEL relates to the 
first of the two families input; i.e., X = f 1 (U) for various V i} and 
therefore determines the positions of the numeric V^ labels. For 
VLABEL = 0, the members of this family are labeled at the end of the 
curve with the greatest value of U. For VLABEL = 1, the members are 
labeled at the other end of the curve, viz, the end corresponding to 
the minimum value of U. For VLABEL - 2, both ends are labeled. The 
numeric labels are automatically offset from the ends of each member. 
The variable ULABEL is similarly defined for the family X - f (V) for 
constant Uj that is either input following the first family or is 
determined by interpolation. However, in order to distinguish between 
the numeric labels for the two families, the position of these labels 
with respect to the end of each curve is determined by the two param- 
eters, XTICK and YTICK, described in the next paragraph. The default 
values for ULABEL and VLABEL are zeroes. 
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The parameters XTICK and YTICK are associated with the var- 
iable ULABEL discussed in the previous paragraph. The numeric labels 
for the members of the family X = fi 2 00 for constant Uj are offset from 
the appropriate end of the curve (defined by ULABEL) by a tick mark the 
length and orientation of which are defined by XTICK and YTICK. An 
example of this may be seen in Figure 1. XTICK is the length of this 
tick mark projected in the abscissa direction and is measured in rasters. 
The raster coordinate system is the lowest level coordinate system 
available on the SD 4060. The abscissa axis is approximately 3100 ras- 
ters long, and a representative value for XTICK is 100, which number is 
used for the default value of XTICK. YTICK is similarly defined for 
the ordinate direction. The ordinate height of the grid is 2300 ras- 
ters, and a default value of 100 is included for YTICK. The numeric 
label for the family X = ^ (V) is automatically set a small distance 
from the end of the tick mark. For all values of ULABEL, a positive 
value of XTICK results in a positive abscissa increment for the tick 
mark. In the ordinate direction, however, a positive YTICK results in 
a positive increment for ULABEL = 0 and a negative increment for 
ULABEL = 1. For ULABEL = 2, a positive YTICK results in a positive 
increment at one end and a negative increment at the other. Tick marks 
may be supressed by setting XTICK = YTICK - 0. If a single orthogonal 
family is read in, there is no option concerning which family is labeled 
with tick marks. However, if both families are input, the second input 
family is labeled with tick marks, and hence, the order of the input 
families may be revised to alter the labeling. 

The three variables DCHAR, NCHAR, and CHFMT define the labeling 
of the vertical scale markers on the plot. The actual scales are de- 
fined by the two variables DELTAU and DELTAV (Section 4.2.1). The scale 
markers are two I-beam symbols and can be seen in the bottom left-hand 
corner of the carpet plot example shown in Figure 1. The horizontal 
positioning of these markers from the left-hand grid margin is one cycle 
for logarithmic abscissa and DELTX for linear abscissa. The vertical 
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positioning is one horizontal grid-line increment from the lowest grid 
margin with a single increment gap between the two I beams. DCHAR is 
a two-dimensional array containing alphanumeric titles for the scale 
markers. DCHAR fl, 1) and DCHAR (2,1) contain the label for the DELTAV 
scale, and DCHAR(1,2) and DCHAR(2,2) contain those for the DELTAU scale. 
Both titles have a maximum length of eight characters. The DCHAR array 
defaults to blanks. NCHAR is the larger of the two numbers that are 
the number of characters in the DELTAU and DELTAV labels. The shorter 
label is filled with sufficient blanks to increase its length to that 
of the longer label. NCHAR defaults to 1. CHFMT is the F format for 
the values of DELTAU and DELTAV written with the scale markers and 
labels. The form of CHFMT is W.D where W is the field width and D is 
the number of decimal places. W is restricted to a maximum of 6, and 
CHFMT defaults to 4.2. 

The final variable XFMT is required for linear abscissa only 
(ABSCA = 0) and specifies the format for the numeric labeling of the 
abscissa scale. The form is W.D where W is the field width and D is 
the number of decimal places. For positive W the format is F, and for 
negative W it is the E format. 
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5. CARPET PLOT GENERATION 

Despite the large number of variables associated with Program CRP, 
generation of a usable carpet plot is not a difficult matter. There are 
three aspects of the finished plot: (1) the matrix of family members, 

(2) the interpolation of the input family to obtain the orthogonal 
family if both are not input and, (3). labeling and titling. It is 
strongly recommended that the first item be treated separately from 
the remaining two items. Further, if interpolation problems occur, 
these should be solved before final labeling and titling are added to 
the plot. 

Obtaining a matrix of family members arranged in such a manner that 
interpolation can be performed is the most difficult part of carpet plot 
generation. Throughout this section the term "usable carpet plot" is 
used. A usable carpet plot is one on which interpolation between mem- 
bers can be performed satisfactorily. A necessary condition for this 
is that members of a given family must not cross each other and that 
members of opposing families must cross only once. In addition, the 
closer each member is inclined at 45° to the horizontal and 90° to 
each member of the orthogonal family, the more easily interpolation on 
the matrix may be performed. It should be stated that there are var- 
iables X = f(U,V) that are not suitable for carpet plot display, but 
many of these can be adequately plotted by restricting the range of U 
and V on a given plot and mapping the entire field with several carpet 
plots. The functions that should be treated in this fashion are those 
for which either X exhibits several maxima and minima or 9X/3U or 
' 9X/3V changes very rapidly in relation to the range of U and V covered. 

The first computer run should be intended to test only the matrix of 
family members. All of the available data should be read in, but only a 
small number of members should be plotted. Initially, only the carpet 
plot control parameters (Section 4.2.1), some of the grid parameters 
(Section 4.2.3), and a small number of labeling and titling parameters 
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(Section 4.2.4) should be specified. If the required plot is logarithmic, 
only the following NAMELIST parameters should be specified: RTYPE, DELTAU, 

DELTAV, UPLOT, NUPLOT, VPLOT, NVPLOT, UMAX, UMIN, VMAX, VMIN, TITLE1, Nl, 
TITLEX, NX, DCHAR, NCHAR, and CHFMT. If the plot is linear, the following 
additional parameters should be specified: ABSCA, XO, XF, DELTX, and 

SUBDIV. If only a single family is input, the parameters RTYPE, UMAX, 

UMIN, VMAX, and VMIN may be removed from this list. All of the other 
NAMELIST parameters discussed in Section 4.2 have their default values 
initially. 


The choice of the six parameters DELTAU, DELTAV, UPLOT, NUPLOT, 

VPLOT, and NVPLOT is important for the initial program run. As noted 
in the previous paragraph, all of the available data should be input. 
However, the parameters UPLOT, NUPLOT, VPLOT, and NVPLOT should be 
used to plot only a small number (six or seven) of the salient members 
of the input family or families. The VPLOT array should include the first 
and last members of the input family X = f 1 (U) for various V^. If the 
orthogonal family is read in, UPLOT should contain the first and last 
members of this family. If the orthogonal family is obtained by 
interpolation, UPLOT should contain values close to the minimum and 
maximum values of U input with the first family X = f (U) . 


Initially DELTAV should be a positive quantity. DELTAU and DELTAV 
not only determine the relative positions of the various family members, 
but also influence the number of horizontal grid lines drawn. For ease 
of interpolation on the finished carpet plot matrix, approximately 40-50 
horizontal grid lines should be drawn. The actual number drawn is output 
on the line printer (Section 7), and the suitability of this number may be 
assessed without the necessity of submitting the SD 4060 plot tape. From 
equations (3), (7), and (9), for the initial computer run, the number of 
horizontal grid lines drawn is, in FORTRAN notation. 


m 


V 


n 


V. 


N = 5 + IFIX (4 + 


DELTAU 


DELTAV 


(ID 



where the subscripts 1 and n refer to the minimum and maximum matrix 
points, respectively. For RTYPE = 1, U m , U : , V n , and V x are input via 
NAMELIST. For RTYPE = 0, assuming that the VP LOT and UPLOT arrays are 
monotonically increasing, V = VPLOT (1) 

and V n = VP LOT (NVPLOT) 

for positive DELTAV. For the purposes of estimating the number of 
horizontal grid lines drawn, it may be assumed that U m is the greatest 
number in the UPLOT array that is less than or equal to the maximum 
U value input for the member X = f L (U) corresponding to V^ = VPLOT (NVPLOT) . 
For example, in Figure 6, VPLOT (NVPLOT) = 18000 n.m., and the maximum 
UPLOT value contained in the Energy array for V^ = 18000 n.m. is E = 4 
MeV, and thus U m = 4.0. In many cases, however, such as that shown in 
Figure 1, U m = UPLOT (NUP LOT) , as each input member contains values for 
the entire range of U. The variable U x may be taken to be the smallest 
number in the UPLOT array that is greater than or equal to the minimum 
U value input for the X = f a (U) corresponding to V^ = VPLOT(l). U t 
will usually reduce to Uj = UPLOT (1) as in Figures 1 and 6. U x and 
U m should be obtained using these definitions by inspection of the 
input data deck. 

Equation (11) is clearly not sufficient to estimate initial DELTAU 
and DELTAV values even if the points (1,1) and (n,m) (Figure 3) and the 
required number of horizontal grid lines are known before the initial 
plot is made. However, first estimates of DELTAU and DELTAV may be 
obtained by assuming that 

U m - ^ V n - V i (12) 

DELTAU = DELTAV 

Assuming that approximately 45 horizontal grid lines are required, 
equations (11) and (12) become 

DELTAU ~ | Um - Uj / 18 (13) 

and DELTAV r | VPLOT (NVPLOT) - VPLOT (1) | / 18 (14) 
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Applying equations (13) and CL4) to the data plotted in Figure 6, 
initial estimates of DELTAU s AE = 0.222 MeV and DELTAV = AALT = 777.77 
n.m. are obtained. It should be reiterated that these estimates are 
obtained before any plots are made. Equations (13) and (14) fix the 
maximum and minimum matrix points on horizontal grid lines. The con- 
ditions that any other matrix point (i,j) will lie on a horizontal 
grid line, which is desirable for ease of interpolation, are that Uj - 
U : and - V x must be integer multiples of DELTAU and DELTAV, respec- 
tively. It is clear that this is not likely to be the case for AE = 
0.222 MeV or AALT = 777.77 n.m., and in general the initial estimates 
of DELTAU and DELTAV should be 'rounded' to more reasonable values. 
Therefore, in the above case, AE is rounded to 0.25 MeV and AALT to 
750 n.m., and the results of the initial computer run are shown in 
Figure 7. It can be seen that sufficient labeling has been added to 
easily identify the family members and scale parameters. A second 
example is shown in Figure 8. In this case, initial estimates of the 
scale parameters were obtained from equations (13) and (14) as DELTAU 
= AB = 0.01333 and DELTAV = AL = 0.0666. These were rounded to AB = 
0.015 and AL = 0.06. Note that in the case of the data shown in 
Figure 8, it was assumed that U x - UP LOT (1) =0.04. As can be seen 
from Figure 8, the minimum matrix point does not, in fact, occur for 
U = B = 0.04. Nevertheless, a usable carpet matrix resulted from the 
use of equations (13) and (14), and this is the initial objective. 

From the two examples shown in Figures 7 and 8, it is evident that 
usable carpet plots can be obtained for the input data. In some in- 
stances this will not be the case, and the various family members may 
cross over each other in some regions of the matrix. This does not 
necessarily mean that the data cannot be carpet plotted. For instance, 
the data plotted in Figure 8 are replotted in Figure 9 with different 
values of the scale parameters, and many of the members cross over each 
other. Attempts at carpet plotting should not be abandoned merely 
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because members cross initially. Two courses of action are possible. 
First, the data may be split into subsets each of which may be 
successfully plotted on separate frames. Second, the scale parameters 
may be varied in order to generate a usable matrix (as in Figures 8 
and 9) . 


Varying the scale parameters should be performed with care since 
they affect not only the ordinate distance for any member but also the 
number of horizontal grid lines drawn. The change in the number of 
horizontal grid lines drawn is, from equations (11) and (13), 


6N ~ - IS 


6 (DELTAU) 6 (DELTAV) 
DELTAU + DELTAV 


(15) 


where the 6 symbol indicates small changes in the independent variable 
following. Further, the change in the vertical extent of a family 
member expressed as the number of horizontal line increments that it 
covers, Nj, is 

^1 „ 6 (DELTAU) Qr _ 5 ( DELTAV) r 

~ w7" DELTAU DELTAV 


for the families X = f : (U) and X = f 2 (V), respectively. Expressed in 
terms of the SD4060 coordinate system, Y s , the change in the vertical 
extent of a given member is 


6Y S „ . 

6 (DELTAU) + ® "1 nT . . 

"6 (DELTAV) v 6N 

Y s 

DELTAU N-lJ 

DELTAV N-l 


(17) 


where 6N is given by equation (15). For positive DELTAU and DELTAV, 
equation (15) indicates that an increase in either scale parameter 
results in a reduction of the number of horizontal grid lines drawn. 
Equation (16) shows that the vertical extent of a given line in terms 
of grid- line increments decreases for a positive change in the scale 
parameters. However, equation (17) for the subject coordinate distance 
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(which is the important variable in the context of modifying the carpet 
matrix) is ambiguous as far as the sign of 6Y S is concerned, since the 
two terms have opposing signs. In general, however, the first term is 
dominant, and for a positive change in a scale parameter, 6Y S is nega- 
tive. For instance, if DELTAU is increased by 10 percent (reducing the 
number of grid lines drawn by one or two depending upon rounding), the 
vertical extent of a given member is reduced by 10 percent in terms of 
grid-line increments but only 5 percent in terms of subject coordinates. 

In Figure 9, the input family is that labeled ' L Values', i.e., 

X = f(B) for various L. The orthogonal family, labeled 'B Values', is 
obtained by interpolation. The members of the input family cross over 
each other primarily because the vertical extent of the higher L-value 
members is too large. From equations (15) and (17), assuming that 
approximately 45 horizontal grid lines are drawn, 5Y S /Y S ~ - 0.5 
6 (DELTAU)/ (DELTAU) . Thus, in Figure 9, in order to reduce the vertical 
extent of the L-value members by half, the value of DELTAU (= AB) should 
be approximately doubled. It should be noted that equations (15) through 
(17) are only approximate, are applicable for small changes in the 
scale parameters, and assume that the maximum and minimum grid points 
do not change. The above expression for <5Y S /Y S is therefore approximate 
and should be used to determine trends only. From equation (15), if 
DELTAU is doubled and DELTAV is held constant, approximately 18 fewer 
horizontal grid lines will be drawn. As a result, DELTAV should be 
reduced to compensate for this. On the bases of the carpet matrix 
shown in Figure 9 and the trends indicated by equations (15) and (17), 
new values of the scale parameters may be estimated as DELTAU = AB = 

0.02 and DELTAV = AL = 0.05. The carpet matrix for these scale param- 
eters is shown in Figure 10, and it can be seen that this matrix is 
acceptable. 

Having obtained a usable matrix, a plot should be made that contains 
all the required members in the matrix. At this point, any interpolation 
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problems should be solved. The two available interpolation schemes, 
SPLINE and CRVFIT, are contrasted in Figures 11 through 13. In each of 
these figures the L-value family is input. In Figures 11 and 12 the 
only members input are those that are plotted. It may be seen that the 
interpolated family obtained with the CRVFIT routine is acceptable 
whereas that obtained with the SPLINE routine has several unacceptable 
variations in the members, particularly for B = 0.04. In general, this 
type of oscillation is likely to occur when a large change of slope 
occurs over a small region of the varying independent variable. If 
interpolation problems are not solved simply by changing the routine 
used, additional input members must be added in order to damp the 
oscillation. For instance, in Figure 12 an input curve for L = 3.2 
may be added. It is not necessary to plot this input curve, and thus 
values of the dependent variable may be input only for those values of 
the independent variable that are contained in the UP LOT array. That 
is, in Figure 12, the L = 3.2 member requires X values for B = 0.05, 

0.1, 0.2, 0.4, and 0.5. The values of X may be estimated from the 
existing plot (Figure 12). The results of this procedure are shown in 
Figure 13, in which it can be seen that the oscillation has been damped. 
The results, however, are still clearly inferior to those obtained with 
CRVFIT routine, and in this case this routine should be used. 

The final phase of the carpet plot generation, the addition of 
titling and final labeling, is a straightforward application of the 
NAMELIST parameters given in Section 4.2.4. As an example, the final 
version of the test plot shown in Figure 7 may be seen in Figure 6. 
Additional family members have been plotted, and the member numeric 
labels for the input family have been moved from one end to the other 
to avoid overwriting some members. In addition, comprehensive titling 
has been added. 
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6. JCL SETUP 

Program CRP is available from the GSFC Program Library in source 
(FORTRAN IV) and object deck (FORTRAN IV H, OPT = 2 compilation) form 
suitable for operation on the IBM 360 series machines. Versions of 
the program for operation on IBM 7094 and UNIVAC 1108 machines are 
also available from the National Space Science Data Center. Because 
of core limitations with this computer, the IBM 7094 version is over- 
layed but offers the same capabilities as the IBM 360 version. 

On the 360/75 machine, compilation of all subroutines requires less 
than 1.0 minute cpu time and 1.5 minutes I/O time. Execution time to 
generate 10 carpet plots with single family input is approximately 1.0 
minute cpu time and 1.0 minute I/O time. Less execution time is re- 
quired if both families are input. The program requires less than 170K 
for the EXEC LINKGO step. 

The JCL required for Program CRP is relatively straightforward and 
is shown in Figures 14 and 15 for source and object decks, respectively. 
The output tape for the SD 4060 plotter is written in 7-track binary at 
556 bpi and may either be user supplied or TEMPSV. The data deck setup 
has been discussed in Section 4 and is shown in Figures 4 and 5. 


Preceding page blank 
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7. PROGRAM OUTPUT AND ERROR MESSAGES 
7.1 Program Output 

A binary, 7-track, 556 bpi SD 4060 plot tape is output by Program 
CRP, The number of frames on this tape is output by the SD 4060 pack- 
age on the line printer. For the case KOLOR =1, this number should 
be 2 plus the number of plots input on the first data card. For KOLOR = 

2, the number should be 2 plus twice the number of plots. For ease of 
publication, it is suggested that 35-mm positive film output be requested 
from the SD 4060. 

Line printer output is provided to enable the user to assess the 
success of the program run before submittal of the SD 4060 plot tape. 

The plotted members of each orthogonal family are output, and a sample 
list is shown in Figure 16. The first family output corresponds to the 
first family read in to the program (X = f(U) for various V^) . The 
second family output corresponds either to the second input family or 
to the interpolated family, depending upon RTYPE. Note that only the 
family members to be plotted are output on the line printer. The 
values of the constant independent variable associated with each member 
input are listed before the plotted members (values of 999.0 in this 
list should be disregarded as these are dummy values). In Figure 16, 
the alphanumeric header for each member and the alphanumeric header 
for the columns of the varying independent variable are taken from 
the DCHAR array input via the NAMELIST. Pairs of dependent and independent 
variables are given on each line-printer line, together with the SD 4060 
subject coordinate of each point (labeled Y) . 

Three checks may be performed on these tables before the plot tape 
is submitted. First, the constant independent variable heading each 
member may be compared to the UPL0T and VPLOT arrays to ensure that the 
intended lines are plotted. Second, the values of Y may be scanned to 
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make sure that all values of Y are in the range 0 to 100 (except the 
999.0 default values), and that the variation of Y for a given member 
is significant. Values of Y outside the range 0 to 100 may result from 
•overhang' of the input curves. That is, the plot is automatically 
scaled in a vertical direction by using the maximum and minimum grid 
points. The input curves may cover a wider range of independent var- 
iables than those encompassed by these grid points, with the result that 
values of Y s > Y(n,m) and Y s < Y(l,l) (see Figure 3) may result. In 
extreme cases, Y s > Y or Y < 0. This problem may be rectified by 
increasing^ the value of PO to allow more space between the bounding 
horizontal grid lines and the bounding grid points in the carpet matrix. 
If the variation of Y for a given member is small (< 10), the result 
will be an almost horizontal line in the plot, making interpolation on 
the matrix an inaccurate process. This may be rectified by making the 
appropriate vertical scale more sensitive, that is, by reducing either 
DELTAU or DELTAV. This should be done with caution, however, since all 
of the members of each family will be affected by this change, and the 
overall effect may not be to improve the plot. The third check that 
may be performed is to examine grid points in both family lists. The 
subject coordinate Y should be identical for a given grid point in both 
lists . 

In addition to the plotted families, the number of horizontal grid 
lines is output. As noted in Section 5, this should be approximately 45 
to facilitate interpolation on the carpet matrix. This number may be 
changed by varying DELTAU and DELTAV as discussed in Section 5. For 
the case of linear abscissae (ABSCA = 0.), the number of vertical grid 
lines is output. 

Under control of the NAMELIST variable IPR (Appendix A), the entire 
NAMELIST may be written on the line printer. For NUPLOT or NVPLOT < 

20, the UPLOT or VPL0T arrays are filled with dummy values of 999.0, and 
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these should be disregarded. Further, when the default values for the 
alphanumeric NAMELIST variables are used, numeric representations of 
these defaults are output on the line printer. These, too, should be 
disregarded. 

7,2 Error Messages 

A number of SD 4060 error codes may be generated by misuse of Pro- 
gram CRP. These are output by the IGS routines associated with the 
SD 4060 plot package. The reader is referred to the SD 4060 manual for 
a complete list of these. However, the two most likely SD 4060 error 
codes are 32 and 33 and refer, respectively, to bad abscissae and bad 
ordinate values. Bad abscissa coordinates are not likely to occur for 
a logarithmic carpet plot for which the abscissa scaling is performed 
automatically. However, if XO and XF are specified incorrectly for a 
linear carpet plot code 32, errors will result. Code 33 errors may 
result from the ’overhang’ of the input curves that has been discussed 
in Section 7.1. These may be corrected by increasing P0 or by ex- 
tending the range of the parameter UPLOT. 

A number of error codes are output by Program CRP and these are 
listed in Appendix C. Section C.l shows four error messages associ- 
ated with NAMELIST input errors, the last of which is a warning state- 
ment only. Program termination is associated with the first three error 
messages. The first message indicates that an incorrect value of RTYPE 
or ABSCA has been entered, and the second indicates that either DELTAD 
or DELTAV has been entered with the unacceptable value of zero. The 
third message indicates a problem with one or more of the five listed 
variables that are required for a linear carpet plot only. This error 
message will result either if these variables have not been specified 
or if XF < XO. The fourth message is a warning indicating that any one 
of the three arrays (TITLE1, 2, and 3) has not been specified. 
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The error messages associated with the interpolation procedure for 
RTYPE = 0 are shown in Section C.2. The two return codes associated 
with the first error message indicate problems with the input data. 

The return code of 1 generally indicates that an input curve is not 
monotonic with respect to the varying independent variable, e.g. , U of 
X = f 0J) for various V = V^. The parameter INT is the number of incre- 
ments of the interpolated curve between grid points and is Np + 1 in 
equation CIO)- The maximum value of INT is M = 48 since the maximum 
length of the interpolated array is 50, and a minimum of two curves 
must be input. The return code of 2 generally indicates either that 
VPLOT is not monotonic or that VP LOT contains a value that was not 
input. IGP is the array of grid points stored for each value of UPLOT. 
Error message numbers two and three result if the UPLOT array contains 
one or more values that are not included in any of the members of the 
family X = f[U) for various V = V^. This will not normally terminate 
the run [unless the entire UPLOT array is incorrect) as the program 
will automatically eliminate these values. The corrected UPLOT array 
and its length are given by error message number three. 
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APPENDIX A. SUMMARY OF NAMELIST INPUT PARAMETERS 


With the exception of the first parameter, IPR, the NAMELIST input 
parameters are defined in order in the four categories of (1) Carpet 
Plot control parameters, C 2 ) Interpolation parameters, (3) Grid 
parameters, and 0*) Labeling and Titling parameters. Each parameter is 
numbered. An alphabetical list of the parameters is given in Table 2 
and this may be used as a convenient cross reference for the list 
given in this Appendix. Section references are given for a more de- 
tailed description of a given parameter. 

No. Name Type Default Section Description 

1 IPR 1*4 None Print Option Flag. 0 for 

NAMELIST variables not printed 
and 1 for NAMELIST variables 
printed. 

2 RTYPE 1*4 0 4.2.1 Input Data Flag. 0 for a single- 

family input and 1 for both 
families input 

3 PO R*4 4. 4.2.1 Determines vertical positioning 

of carpet matrix. Number of 
horizontal grid increments below 
minimum grid point and above 
maximum grid point. Should be 
whole number. 

4 DELTAU R*4 0. 4.2.1 Increment in U (the independent 

variable not held constant for 
the members of the first family) 
corresponding to one increment 
between successive horizontal 
grid lines 

5 DELTAV R*4 0. 4.2.1 Increment in V corresponding to 

one increment between successive 
horizontal grid lines 
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No. 

Name 

Type 

Default 

Section 

Description 

6 

VP LOT 

R*4 

None 

4.2.1 

The array of Vs for which 
X = f(u) are plotted. VPLOT 
e V^, where Vj^ are the values 
of the independent variable 
held constant for the members 
of the first family input. 

7 

NVPLOT 

1*4 

0 

4.2.1 

Length of the VPLOT array. 
Maximum 20 

8 

UP LOT 

R*4 

None 

4.2.1 

The array of Vs for which 
x = f (V) are plotted. UPLOT 
£ Ui for RTYPE = 1. 

9 

NUPLOT 

1*4 

0 

4.2.1 

Length of the UPLOT array. 
Maximum 20 

10 

CHOICE 

R*8 

'CRVFIT' 

4.2.2 

Interpolation Method Flag. 
Applicable for RTYPE = 0 only. 
'CRVFIT 1 for subroutine CRVFIT 
and 'SPLINE' for subroutine 
SPLN2 and function SPLIN(Y) 

11 

UMAX 

R*4 

0 

4.2.3 

Maximum grid U value for input 
data. Omit for RTYPE = 0. 

12 

UMIN 

R*4 

0 

4.2.3 

Minimum grid U value for input 
data. Omit for RTYPE = 0. 

13 

VMAX 

R*4 

0 

4.2.3 

Maximum grid V value for input 
data. Omit for RTYPE = 0. 

14 

VMIN 

R*4 

0 

4.2.3 

Minimum grid V value for input 
data. Omit for RTYPE = 0. 

15 

KOLOR 

1*4 

1 

4.2.3 

Flag for number of SD4060 frames 
output. 1 for single-frame 
output. 2 for double-frame 
output 

16 

ABSCA 

R*4 

1 

4.2.3 

Flag for type of abscissa scale. 


0 for linear and 1 for logarithmic 
scale 
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No. 

Name 

Type 

Default 

Section 

Description 

17 

XO 

R*4 

None 

4.2.3 

Minimum abscissa value in input 
data units. Omit for ABSCA = 1. 

18 

XF 

R*4 

None 

4.2.3 

Maximum abscissa value in input 
data units. Omit for ABSCA = 1. 

19 

DELTX 

R*4 

None 

4.2.3 

Increment in input data units for 
which vertical grid lines, tick 
marks, and labels are given. 

Omit for ABSCA = 1. 

20 

SUBDIV 

R*4 

1 

4.2.3 

Number of increments into which 
DELTX is subdivided by drawing 
vertical grid lines. Omit 
ABSCA = 1. 

21 

TITLE1 

R*4 

Blanks 

4.2.4 

First line of title in page header 

22 

N1 

1*4 

1 

4.2.4 

Number of characters in TITLE1. 
Maximum 60 

23 

TITLE 2 

R*4 

Blanks 

4.2.4 

Second line of title in page header. 

24 

N2 

1*4 

1 

4.2.4 

Number of characters in TITLE2. 
Maximum 60 

25 

TITLE 3 

R*4 

Blanks 

4.2.4 

Third line of title in page header 

26 

N3 

1*4 

1 

4.2.4 

Number of characters in TITLE3. 
Maximum 60 

27 

XTITLE 

R*4 

Blanks 

4.2.4 

Abscissa title. Appears below 
abscissa scale. 

28 

NX 

1*4 

1 

4.2.4 

Number of characters in XTITLE. 
Maximum 60 

29 

LABE LI 

R*4 

Blanks 

4.2.4 

Alphanumeric label for one family 
of curves 

30 

NLAB1 

1*4 

1 

4.2.4 

Number of characters in LABEL1. 
Maximum 36 

31 

XPOS1 

R*4 

0 . 

4.2.4 

Horizontal position of first 


character of LABEL1 in inches. 
Measured from left-hand bounding 
grid line. -1.25 £ XP0S1 <_ 8.0. 
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No. 

Name Type 

Default 

Section 

Description 

32 

YP0S1 R*4 

0. 

4.2.4 

Vertical position of first 
character of LABE LI in inches. 
Measured from lower bounding 
grid line. -0.4 < YP0S1 £ 5.9. 

33 

THETA1 R*4 

0. 

4^ 

» 

to 

• 

Orientation of LABEL1 in degrees. 
Measured counterclockwise from 
the horizontal. 0. < THETA1 < 90. 

34 

LABEL2 R*4 

Blanks 

4.2.4 

Alphanumeric label for other 
family of curves 

35 

NLAB2 1*4 

1 

4.2.4 

Number of characters is LABEL2. 
Maximum 36 

36 

XP0S2 R*4 

0. 

4.2.4 

As XPOS1 bar for LABEL2 

37 

YPOS2 R*4 

0. 

4.2.4 

As YP0S1 bar for LABEL2 

38 

THETA2 R*4 

0. 

4.2.4 

As THETA1 bar for LABEL2 

39 

ULABEL 1*4' 

0 

4.2.4 

Flag for numerical labels for lines 
of constant llj . 0 for end of 

line with maximum V value. 1 for 
end of line with minimum V value. 

2 for both 

40 

VLABEL 1*4 

0 

4.2.4 

Flag for numerical labels for 
lines of constant V^. 0 for end 

of line with maximum V value. 1 
for end of line with minimum U value 
2 for both 

41 

XTICK R*4 

100. 

4.2.4 

Horizontal distance in rasters for 
tick mark used with numeric labels 
for the family X = f(V) for various 
Uj . Horizontal distance positive 
for positive XTICK. 


38 



Name Type Default Section 


Description 


YTICK 


DCHAR 


NCHAR 


CHFMT 


XFMT 


R*4 100. 4.2.4 Vertical distance in rasters for 

tick mark used with numeric labels 
for the family X = f(V) f° r 
various Uj . Vertical distance 
positive for positive YTICK and 
ULABEL = 0. Vertical distance 
negative for positive YTICK and 
ULABEL = 1. 


R*4 Blanks 4.2.4 


1*4 1 4.2.4 


R*4 4.2 4.2.4 


R*4 2.0 4.2.4 


Array containing alphanumeric 
labels for the scale markers. 

DCHAR (1,1) and DCHAR (2,1) contain 
label for DELTAV. DCHAR (1,2) and 
DCHAR (2, 2) contain label for 
DELTAU. Both labels have a 
maximum length of eight characters. 

The larger of the two numbers 
of characters in the DELTAV and 
DELTAU labels 

Print format for numeric values 
of DELTAU and DELTAV. The form is 
W.D, where W is the field width 
(W <_ 6) and D is the number of 
decimal places. 

Print format for numeric labeling 
of abscissa scale. The form is 
W.D, where W is the field width 
and D is the number of decimal 
places. For positive W the 
format is F, and for negative W 
the format is E . 



APPENDIX B. PROGRAM LISTING 


DIMENSION Y 1 { SO *20 > * Y2 { 50* 20 ) 

Q I MENS ION XP< 50 ) • YP( SO ) , XV (50 ) *Y V< 50 > 

DIMENSION NPTSI20 ) *LGP < 20) 

DIMENSION X0AR <6 > * YBAR < 6 ) 

D I MENS ION XHL ( 20)i YHL' ?0)*XVL(2Q)*YVL(20),PVX{2) *P V Y C 2 > 

DIME NS I ON XVL B ( 20 ) *YVLQ{ 20 > ♦ XHR( 20 > • YHR < 20 ) 

DIMENSION CHAR4<2*2) 

REAL L * LP 

R£AL*a CHOICE 
REAL 40 CHARI* CHAR2 

*_0GICAL*1 EQUAL « BLANK »LAB 1« LA02 

COMMON /VAL/ BP (20 ) • Ql 50 #2 0 ) *F l ( 50 * 20) *E2( 50 ,2 C ) *1 GP <20 ) « T ( 20) • 

1 L ( 50 * 20 ) , LP (20 

integer rtyp.elabel 

COMMON / IN I T/ DELL *D EL E3 * PO * NL * NF P t N I P ♦ ID I M * RTYP 

COMMON /SCAL / CHO ICE* A BSC A * XG *XF * DCLT X* XFMT»RLMAX* 6MA X * RLM IN# 

1 BMIN*YY • KOLOR « NLP* NBP 

COMMON / T I TL/ N 1 *T IT LE I ( 1 5 > * N2 • T l TLE2 C 1 5 > • N3 , T I TLE 3 i 1 5 ) * NX ♦ 

1 XTITLE{15>*NHT.H0RIZ(9)*NVT,VERT(9)*XHZ»YH2,THETA1» 

2 X VT * YV T, THET A 2 *DCH AR (2*2) t NCH AR* CHFMT * BLABEL ,LLA DEL , 

j DELTA. DELTA! *SURDI V 

COMMON /PLOT/ A <200) 

COMMON /PSI/ PSl l < SO, 20 ) * PS 12(50 , 20 > 

COMMON /L A8G/ L A B l ( 8 ) • L AB2 ( B ) 

EQUIVALENCE <P5 I2( 1* l ) « Y2< 1 * 1 W • ( PSU ( 1 ♦ l » » Y 1 U , U ) 

EQUIVALENCE (GEAR l ,L AD 1 ( l > ) , I CHAR2 , L A0 2 ( l ) ) 

EQUIVALENCE (CFAR4( 1 . 1 ) , CH AR 1) , ( CH AR 4 ( 1 ♦ 2 ) ,CHAR2 ) 

DATA EQUAL/1 H^/, BLANK/ 1H / 

INDEX! I) = ( I— 1)/I INI + 1 

NCURVE-20 

IDIM = 50 

RE AO ( 5* 1007) KARP 
CALL MODESGC A*0 ) 

START NErf PLOT 


l CALL I N IT AL 
CALL READ 
II NT = IDIM/NL 
IF (PTYP.LT.2) GO TO 10 
NHL - NLP 
NVL = NOP 
GO TO 40 

10 IF C KT YP • E Q • 1 ) GO TO 20 
NIP = NBP 
NFP - NLP 

CALL FtNDF(BP.E,Fl,F2#T) 
CALL L INT (L.T*F2) 

CALL GR I DP T (LP»L* I GP ) 

NHL = NFP 
NVL = NIP 
GO TO 40 
20 NIP=NLP 
NFP^NBP 

CALL F INDF (LP.L. F2*F l, T ) 
CALL L I NT C D*T*F 1) 

CALL GKIDPTI BP, B. I GP ) 

NHL - NIP 
NVL = NFP 
40 K = 1 

DO 1 00 1 “I , NHL 

CONST = LPC I ) /DELL 


/ 


Preceding page blank 
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IF ( RT YP#GT ♦ 0 ) GO TO 45 
N= INDEX! IGP! I ) > 

K ■= N 
GO TO 50 
45 N = I 

IF {RTYP # EQ • 2 1 K = N 
50 DO 100 J=1 * IDI* 

P = 0 ( J * K ) 

IF i P*EQ*999 • *GR*F 1( Jt N) «6G« 990* ) GO TC 100 
PS I l { J * N ) “ CONST ♦ P/DELB 

100 CONTINUE 
K = 1 

DO i 05 1 = 1* N VL 

CONST = BP < I 1 /DEL 8 

IF !RTYP*NE#ll N=I 

IF ( RT YP«EQ» 1 ) N= 1 NDEX ( I GP C 1 > ) 

IF (RTYP*EG*l> X=N 
IF !RTYP.EQ*2> I 
DO 1 05 J=1 * ID1W 
P = L( J#K> 

IF ( Pm EQ *999 * *0R *F2( J* N ) *E 999* > GO TC 105 
PSI2!J*N) = CONST ♦ P/OELL 
105 CONTINUE 

IF ( HT YP«LT • 2 ) GO TO 110 
PMX - HLMAX/DELL + EMAX/DELD 
PM N = RLMIN/DEL L + BMIN/DEL9 
RNUM= I NT ! 2 • *PO + PMX - PMN) 

PHI = PO - RLM IN/DELL - BM IN/D EL B 
GO TO 115 

110 LF { RT YP * EG • 0 ) CALL SCALE! P5I2 «RNU^*PHI) 

IF ! RT YP • EG* 1 ) CALL 5CAL6(PSI1* RNUP#FH1) 

115 CONTINUE 

I RNUM= RNUM + l • 

IrfR IT EC 6*10 17) IRNUM 

IF CAOSCA * EQ * 1«) GO TO 116 

NQLlN£=< ( ( XF-XO)/DELTX >*SOBO IV** 00 01 ) + 1 

HRITEC6. 1003) NOLINE 

116 CONTINUE 

DIVN - YY/RNUM 

DO 120 I L= 1 * NHL 

IF ! RT YP • GT # 0 ) LN = IL 

IF (RTYP.EQ40) LN-INDEX! IGPC IL I) 

DO 120 10=1* IDEM 
PI = PSI 1 ( IB.LN) 

IF ( PI • NE* 999 • ) Y 1 ( I B • LN > = D I VN * ( PH I + P 1 > 

120 CONTINUE 

DO 135 IB=t*NVL 

IF (RTYP*E0.1> LN= INDEX! tGPC 10 )) 

IF («TYP*NE*1 ) LN = IB 
DO 135 IL=l*lDIM 
P2 s PSI2! IL • LN ) 

IF CP2.NE. 999.1 Y2(IL*LN) = D I VN * < PH I ♦ P2 ) 
135 CONTINUE 

IF ( RT YP .EQ. 2 ) GO TO 129 
I ERR = l 
I ERR » 0 
1= 1 

INVRT1 =0 
I N VRT2=0 

IF ( RT YP • FQ • 11 GO TO 137 0 

1360 IF ! Y2 (1*1) *N6 ♦ 999.) GO TO 136 1 

1=1 + 1 

GO TO 1360 

1361 IF !Y2( 1 + 1*1 ) *GT * Y2(l*l)l GO TO 1330 
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CALL REVERS(F2*Y2,LGP) 

INVRT1=1 

DO l 362 3=1.20 

IGP( J >=LGP< J > 

1362 CONTINUE 

GO TO 1300 

1370 IF { Yl (1*1) .NE « 9 99 • ) GO TO 1371 
1 = 1 + 1 

GO TO 1370 

1371 IF (Y1(I+1,1) .GT • Yl(I.U) GO TO 1380 
CALL REV ERS< F 1 » Y l • LGP ) 

INVHT2=1 
DO 1372 J= 1 » 20 

IGP{ J)=LGF( J) 

1372 CONTINUE 
1 380 CONTINUE 

IF ( RT YP.EQ.0> CALL I NT ERP < Y 2 • F2 ,IGP,IINT,50.20.IERR> 

IF (RTYP • EQ * 1> CALL I NTE RP ( Y l * F 1 • l GM * X I N T , 50 • 2 0 * IERR > 

IF (1NVRTI »EQ * 0) GO TO 1290 
CALL REVFRS<F2.Y2,LGP) 

DU 1385 J= 1*20 
' IGP( J)=LGP< J) 

1385 CONTINUE 

GU TO 1399 

1390 IF ( 1NVPT2 .ECU 0) GO TO 1399 
CALL REVERSCF l *Y 1 ,LGP) 

DO 1395 J = 1*20 

IGP( J)=LGP( J > 

1395 CONTINUE 

1399 CONTINUE 

137 IF ( IERK «£Q«0) GO TO 129 
WRITE (6*20011 IERR 

129 CHAR4C I.1)=DCHAR( 1*1 ) 

CHAR 4 ( 2* U=DCHAR( 2, 1 ) 

CHAR4< 1* 2)=DCHAR< 1 .2 J 
CHAR4( 2*2) =DCH AR (2*2) 

DO 130 1=1*8 

IF (LA81(I ).NE. EQUAL) GO TO 130 
LABI ( I ) = BLANK 
GO TO 132 

130 CONTINUE 

132 DO 133 1=1*0 

IF ( LA B2( I ) . NE .EQUAL 1 GO TO 133 
LA 62 ( i > = BLANK 
GO TO 140 

133 CONTINUE 

140 WRITE! 6* 1000 > CHARI 
IN = 1 

DO 150 1=1, NHL 

WRITE! 6* 1001 > CHAHULPt I ) * CH AR 2 * CH AR 2 
IF ( RT YP *E Q • 0 ) INO= I NO EX ( I GP ( 1 ) ) 

IF t RT YP *GT«0 ) IND = I 
IF (RTYP.NC.l) IN = SND 
DO 150 J=1»1DIM*2 
K = J+ 1 

IF (FI < J* I ND) .NE .9 99 ..OR.F1 < K, IND) .NE .9 99. ) 

*W«1TE( 6* 1002) (FKN, IND), Y1(N, I NO >,d(N, IN ) *N= J*K) 

150 CONTINUE 

WR ITE< 6* 1000 > CHAR 2. 

IN = l 

00 200 1=1, NVL 

WRITE! 6, 1001 ) CHAR 2* HP < I U CHAR I. CHARI 
IF ( RTYP.NF,. 1 ) IN 0 = I 
IF (RTYP.EG.l) I N D = l NO E X ( I GP { I > ) 

IF (RTYP.NE.O) IN = IND 
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DO 200 J=1 f ID IM*2 
K = J+ 1 

IF (F2( J, I ND > ♦NE.999 * .QR.F2( K* I NO > . NE . 99 9 . ) 
tMR X T E< 6*1002) <F2(N. IND) . Y2(N. INO) .LtN.U) . N=J.K) 

200 CONTINUE 

** BEGIN PLOTTING 

I PASS = 0 

220 CALL RSETMG(A) 

IP ASS = IPASS 4 1 

IF < IP ASS. EQ » 2 ) CALL P AGEG (A*0*1*1) 

CALL S£TSMG< A* 19.4095. > 

CALL SET SMG( A. £0*3071. ) 

RMIN = 500. 

RM AX = 3600. 

YMIN ^ 2 00. 

YM AX =2500. 

XI NCH=8.0625 
Y I NCH= 5 . 90 

X«PERI=(RMAX-RMIN>/X INCH 
YRPERI~t YMAX-YHN )/Y INCH 
XVT1=XVT*XRPER I+RMIN 
XHZ i=X HZ 4XRPER I + RM IN 
Y VT l=Y VT^YRPER I+YM IN 
YHZ 1 =YHZ*YRPE R I4YM IN 

IF (K0L0R.EQ.2) CALL SUPER <RMlN* RMAX.YMIN *YMAX ) 
CALL OBJCT G< A* RM IN*YMI N* RM AX . YMAX) 

IF < IPASS . EQ . 2 ) GO TO 200 
XMAX = FUM) 

XMIN = XMAX 
DO 250 1=1 .NHL 

IF (RTYP.EQ.O) IND = I NDEXC I GP ( 1 )) 

IF (RTYP.GT.O) IND = I 
NPTS< I ) = 0 
DO 240 J= 1 , IDIM 
X = FI 1 J* IND ) 

IF (X.EQ.999.) GO TO 240 
IF (X.LT.XMIN) XMIN = X 
IF {X.GT.XMAXI XMAX - X 
NPTS< I ) = NPTS( 1)41 
240 CONTINUE 
250 CONTINUE 

QO 260 J=1.NCURVE 

DO 255 1=1* IDIM 



IF 

(F2( I . J> 

. EQ . 

999. ) 

GO TO 255 


IF 

<F2C I* J) 

• GT. 

XMAX ) 

XMAX=F2( I* J) 


IF 

(F2( I, J) 

• LT. 

XMIN) 

XM i N=F 2 ( I* J) 

255 

CONTINUE 





260 CONTINUE 

300 IF ( A8SCA. EQ. 1 . > GO TU 400 
CALL SUGJEGC A* XO . 0 .* XF , YY ) 

CALL SET SM G(A*23.0#} 

IF < IPASS .NE* 1) GO TO 350 

D 1 VS N= DELTX/SUEDIV 

CALL GR IOG( A.D IVSN ,D IVN,0. 0) 

350 IF ( KOLOR .EQ.2.AND.I PAS5*EQ. 1 ) GO TO 220 
XS I Z E = — 30 . 

XPL ACE-=— 1*0 

CALL SETSMGC A* 102*XSIZE) 

CALL L A6ELGC A • 0 * OELT X* 0 * XFMT ) 

CFAC = (RMAX-RMINI/I XF-XO) 

GO TO 455 
400 CONTINUE 
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IF ( IPASS.EQ.2) GO TO 430 
IF ( XM I N«LT# l . ) XMIN-1* 

I 1 = ALQG I 0 ( XMIN > 

I 2= AL G G 1 0 < XMAX 1*1 
XO =10 .**1 i 
XF = 10 .**12 

CFAC = ( RMAX-RMIN )/( ALGG10<XF >-ALQGl 0< XC) > 

430 CALL SUOJEGC A, *0,0 ., XF ,YY> 

CALL SETSMG< A, 22, 1.) 

CAlL SETSMGI A. 100. 2. ) 

0=2. 

C = DIVN 

IF ( [PASS. EO* 1 ) CALL GR I OG ( A , D ,C , 1 , 0 > 

IF ( KOLOR.EQ.2 .AND. I PASS.EQ. 1 ) GO TO 220 
YLAB = -YY/60 . 

NC2=2 

I 1 0=10 

DO 450 IX=1 1 * 12 
NC = l 

IF { IX .GT* 9 ) NC = 2 
IF ( IX.GT.99) NC = 3 
XLAB = 10* ** I X 

IF (IX.EQ.I21 XLAB = XLAB - NC *5 0 . 4 9 . * 1 0 * * 4 < 1 2- l ) / CF AC 
XI GLAB=XLAB-< 125./CF AC UXLAB 
Y 1 0LAB=YLAB- 1.0 

CALL NUM0RG { A , X 10LAB ,Y 1 OLAG, NC 2, I 10 ) 

450 CALL NUMBRG( A, XLAStYLAE* NC * I X ) 

455 CUNT INUE 

IF (KOLOR.EG.l) GO TO 450 
CALL SET SMG< A , 1 13, RM IN ) 

CALL SETSMGI A, 1 14, YM IN 1 
CALL SETSMGI A, 1 15t RM AX ) 

CALL SETSMGt A , I 16. YMAX ) 

456 CALL SETSMGC A , 45, 1 .5 > 

XPL ACE=-2. 5 

IF ( ABSCA .EQ. 0.) XPLACE=-1.5 
CALL SETSMGt A, 104.XPLACE) 

CALL T ITLEGt A, NX, XT! TLE, 0# DUMMY, 0, DUMMY) 

X PL ACE = 2.0 

CALL SETSMGC A , 104, XPLACE > 

CALL T ITLEGl A, 0, DUMMY, 0, DUMMY, N3 , TI TLE 3 > 

XPLACE = XPLACE + 1.5 
CALL SET SMG ( A, 10 4, XPLACE > 

CALL T ITLEGC A, 0, DUMMY, 0 , DUMMY , N2 , TI TLE 2) 

XPLACE = XPLACE + 1.5 

CALL SETSMGC A, 104, XPLACE ) 

CALL T ITLEG( A»C« DU MM V , 0 , DU MM Y , N1 .TITLE 1 ) 

CALL SETSMGC A, 30, 2.0 ) 

PLOT FIRST FAMILY 

BACK = RM IN - 200. 

YF AC = ( YMAX-YMIN )/YY 

DELX = l .5*0ELTX 

IF (ABSCA .EQ. 1.) A LO GXO= AL OG l 0 ( X 0 ) 

DU 500 L INE= I * NHL 

IF (RTYP.EQ.O) IND=INDEX( IGP(LINE) ) 

IF (RTYP.GT.O) IND = LINE 
LIM = NPTS(L INE) 

DO 460 NP= 1 , 10 IM 
X = FI (NP. INO ) 

Y = YX ( NP, IND ) 

IF < X.EQ.999. .GR.Y .EQ. 999. ) GG TO 460 
XP(I) = X 
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yp( n = y 

NPT = 1 
II = NP ♦ 1 

GO TO 470 
460 CONTINUE 
4-7 0 L I M = LlM + 11-2 
DO 4 75 NP- l 1 « L IM 
NPT = NPT + 1 

XP { NPT ) - FI ( NP* IND) 

YP(NPT) = Y1 ( NP* IND) 

IF < XP< NPT > .GE.XO .AND.YP< NPT) .NE.999.) GO TO 475 
IF { (XP(NPT-l) - XP ( NPT ) ) • GT *0 EL X ) GG TO 474 

YP( NPT ) = YPINPT-1) ♦ (XO - XP(NPT-l)) / (XP<NPT> - XP(NPT-1>> * 
1 ( YPCNPT) - YP ( NPT- 1 )) 

XP(NPT) - XO 
GO TO 475 

474 NPT - NPT- I 
GO TO 476 

475 CONTINUE 

476 CALL L INESG(A»NPT.XP»YP) 

460 XR1GHT=XP( 1 ) 

YR I GHT — YP{ 1 ) 

XL £FT=XP(NPT ) 

YLt£FT=YPCNPT ) 

IF (XRIGHT . GT . XLEFT) GO TO 466 
XR I G HT =XP( NPT I 
YR iGHT = YP( NPT ) 

XLHFT=XP( l > 

YLEFT=YP( 1 ) 

486 IF C ABSCA.EO. I . ) GG TO 490 

XHL(LINE) = BACK + ( XL EFT- XO > *CF AC 
XHR(LINE) = RM IN * ( XR I GHT-XO ) *C F A C + EO. 

GO TO 495 

490 XHL(LINE) = BACK + ( ALOG IOC XLEFT J — ALOG XC ) *CF AC 

XHR ( LINE) - RM IN + C AL OG 1 0 C XR I GH T > - ALO G XC ) 4 CF AC + 10. 

495 YHR(LINE) = YM IN + YRIGHT*YFAC 
500 YHL(LINE) = Y M IN ♦ YLEFT4YFAC 

PLOT SECOND FAMILY 

DO 600 LINV=1*NVL 

DO 550 NP= 1 * IDIM 

IF (RTYP.NE*!) IND = L I NV 

IF { RT Y P • E Q ♦ 1 ) IND= INDEX! I GP (L IN V) ) 

X = F2INP* TND) 

Y - Y2INP* IND ) 

IF { X.EQ.999 • .Ott.Y .EQ. 999. .OR. X.LT. XG. CR -X.GT. XF ) GO TO 550 
KT * NPfl 
XV < 1 J = X 
YV ( 1 > = Y 

NPT = 1 

GO TO 520 
550 CONTINUE 
520 DQ 525 NPL T = K T * IDIM 
X as F2( NPLT# INC) 

Y = Y2 ( NPLT » IND) 

IF IX. EQ. 999. . OR. Y. EO. 999. .OH. X.LT.XO.CR.X.GT.XF) GG TO 530 
NPT = NPT + l 
XV ( NPT I = X 
YV(NPT) = Y 
525 CONTINUE 

530 CALL L INESG( A * NPT » XV . Y V ) 

540 CONTINUE 
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IF t ABSC A* EQ. 1*1 GO TO 542 
XVL(LINV) = RMIN 4 ( XV ( NPT ) -XO >4 CF AC 
XVLBCLINV) - RMIN 4 ( X V ( 1 ) -XO ) *C FAC 

GO TO 545 

542 XVLCLINVI = RMIN + C AL OG 10 C X V C NP T ) ) - ALOGXG ) *CF AC 
XVLBCLINV) = RM IN 4 ( ALOG1 0( XV ( 1 ) ) - ALOGXO>*CFAC 

545 YVLCLINV) = YMIN 4 YV(NPT)*YFAC 
600 YVLBCLiNV) = YMIN 4 YV(1>*YFAC 

LABEL ANO TITLE 

CALL SETSMG (A«23»0.) 

CALL OBJ CT GCA.0..YMIN* 40 95 •? YMAXl 
CALL SUBJEGC A.0.*YMIN.4095.«YMAX) 

CALL SETSMGC A* 14# 1 • > 

IF (ABSCA .EG. 0.1 X 0 AR ( 3 ) = RMIN 4 CF AC *DEL TX 

IF ( ABSCA. EO. 1.1 XEARC3) = RMIN 4 AL QG 1 0 < 2 • ) *CF AC 
XBAR(l) - X0AR(3) - 25. 

XB AR < 2 ) = XBAR13) 4 25. 

XB AR (4 ) = XB AR 13 ) 

XBAR(5) - XBAfiCl) 

XB AR C 6 ) = XBARC2) 

CALL SETSMG ( A* 30*2*0 ) 

UNIT = DIVN4YFAC 
DO 675 J = l*2 
AD * 1 • 

IF (J.EQ.2) AD = 3. 

DO 650 M=l .6 

IF CM.LE.3) Y E AR ( M ) = AD*UNIT4 YMIN 

650 IF < M • G T* 3 ) YBARCM) = CAD4l.)*UNlT + VMIK 
CALL L INESGC A*G*XBAR*Y BAR > 

YCOR = YBAR( 1 ) 4 DIVN/2.4YFAC 

XL - X BARI 2 1 4 40. 

CALL SETSMGC A. 52* 2.) 

CALL LEGNDG( A. XL. YCOR. W 1HV) 

CALL SETSMGC A. 52.0.) 

CALL TEXTGC A.NCHAP *OCHAR( 1 * J > ) 

CALL GET SMGC A. 1 7* XL A ST ) 

CALL GETSMGC A. 18* YLAST ) 

XX * XL AST 4 50. 4 47.+NCHAR 

DL ~ ABS C DELL ) 

OB^ABSC DELB ) 

IF CJ.EO.l) CALL NUM BRG < A « XX * YLA ST » C HP M T *DL ) 

IF (J.EQ.2) CALL NUM BR G ( A * XX * YLA ST • C HF MT » D B 1 
675 CONTINUE 

CALL SETSMGC A*45. 1 -0 ) 


DO 

700 

1 = 1* 

NHL 


P- 

ABSCLPC l ) 1 


IF 

(P 

.LT. 

• 01 ) 

FMT=e. 3 

IF 

< P 

.GE. 

• 0 1) 

FMT=5. 2 

IF 

( P 

.GE. 

10. ) 

FMT-E. 1 

IF 

( P 

.GE. 

100 • 

) FMT-6.1 

IF 

<P 

.GE. 

1000 

. ) FMT=6.0 

IF 

( P 

.GE. 

10000. > FMT=7.0 

IF 

C P 

.EQ. 

o 

• 

o 

FMT=4. 1 


IF C LLABEL #NE » 1 > CALL NUMB RG { A * X HL < I ) * YHL ( I > *F NT *LP ( I ) ) 
700 IF ( LLABEL #N£ . 0 ) CALL NUMB RG < A .X HR < I ) * YHR ( I) «F MT »LP ( I ) ) 
CALL SETSMGC A .46, THETA 1 > 

CALL LEGNDGC A* XH21 tYHZ 1 • NH T » HQ R I Z > 

CALL SETSMGC A*46.0.) 

DO 800 I=1*NVL 
V~ ABSC BP ( I ) ) 

IF C V .LT. .01 ) FMT‘6. 3 
IF (V .GE. .011 FMT=5.2 
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IF 

( 

V 

« GE . 

to 

• ) 

FM T — 

5* 

1 

IF 

! 

V 

*GE. 

10 

0.) 

FMT 

=6 

• 1 

IF 

( 

V 

• GE . 

10 

OC . 

) FM 

T= 

6.0 

IF 

( 

V 

• GE * 

10 

000 

.) F 

MT 

= 7,0 

IF 

( 

V 

* EQ • 

0* 

0 ) 

FMT- 

4. 

1 

IF 

( 

BLABEL « 

• 60 

. n 

GC 

TO 

750 


PV X ( 1 ) = XVL( I) 

PVX(2) = XVL(I) ♦ DELTA 
PV Y t n = YVL( 1 > 

PVY(2) = YVL(I) + DELTA! 

CALL L I NESG( A» 2tPVX» PV Y ) 

DX=40. 

DY = 30 ■ 

X=PVX! 21 -S I GN( CX, DELTA ) 

Y=PVY< 2>+S IGN< DY* DELTA 1 ) 

CALL NUMBRGI A . X*Y,FMT* BP( l ) ) 

IF ( BLABEL •EQ #0 ) GO TO 800 
750 PVX< 1) - XVL6C I ) 

PVX(2) = PVX(I) + CELT A 
P V Y ( l ) = YVLB< I ) 

PV Y ( 2 ) = P VY ( I ) -DELTA1 
CALL LINESG( A*2*PVX,PVY) 

OX-4 0 * 

DY —3 0 * 

X=PVX( 2)-SIGN{CX*DELTA ) 

Y=PV Y ( 2) -SIGN ( CY, DELTA X) 

CALL NUMbRG( A*X#Y*FMT* BP( I >> 

800 CONTINUE 

CALL SETSMG! A * 46* THETA 2 ) 

CALL L£GN0G< A ♦ X VT 1 , Y VT 1 * NV T * VE R T > 

CALL SET SMG< A, 46*0.) 

CALL PAGEG< AiOiltl) 

KARP - KARP - 1 

IF (KARP.GT.O) GO TO 1 
CALL EXITGCA) 

STOP 

1000 FORMAT! 1H1 *55X *10HLI NES OF CONSTANT ,A6 /> 

1001 FORM ATI 5BX*A0 * 1 H- *F 10 • 3/ l 3X » 7HDEP V AR « 9X • l H Y * 1 1 X • A 6 ,40X*7HDEP VA 
*R , 1 G X * 1HY* 1 0 X * A8 /> 

I 00 2 FORMAT ( 1 OX, 1PE 12.5.4X* 0PF6 . 1 *5X* IPE12.S*37X,1PE12.5*4X,0PFG.1*5X* 

♦ 1 PE 12.5) 

1003 FORMAT! 1H0* 3X.21HNC. OF VERT GRID LINES CRAWN - ,12) 

1007 FORMAT! I2> 

1017 F O HM AT ( l HO * 3 X » 32 HNC * OF HORIZ GRID LINES DRAWN = ,13) 

1099 FORMAT(5X, I4,ieF12.2)> 

2001 FORMAT ( 1H1 ,////// * 36X* 

1 6QH * ♦*** AN ERROR CONDITION HAS BEEN DISCOVERED IN INTERP ****** 
2////,49X*26HTHE ERROR RETURN CODE IS *11, 7H* WHERE,//, 

3 44X , 40 HI — INT IS LESS THAN 1 CR GREATER THAN M*/* 

4 4 4X « 44H (THIS MEANS INT CAN NOT BE A VALID INDEX*/, 

S44X*27H INCREMENT IN X I N , Y IN ) . * / / . 

644X*52H2= AN INDEX IN IGP(N) IS OUTSIDE THE INDEX RANGE OF,/* 
744X.50H THE FIRST THROUGH LAST DATA POINTS OF CURVE N.) 

END 

BLOCKDAT A 

INITIALIZE NAM EL I ST ( SEE ALSO INlT) 

INTEGER RT YP * EL ABEL 

COMMON / t N I T/ DELL *DELB *PO *NL * NF P*NIP, IDIMtRTYP 

COMMON /SCAL / CHO ICE, ABSCA* XO ,XF »OELT X * XF M T ,RLM A X « BMA X * RL M 1 N* 

1 BM IN # YY • KOLOR , NLP • NBP 
REAL *8 CHOICE 
REAL L * LP 
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o> ir Ui w upo 


COMMON /VAL/ EP < 20 I • BI 5 0 « 2 0 ) • F 1 < SO . 2 0 ) « F 2 < 5 0 *2 0 ) • I GP ( 2 0 ) * T ( 2 0 ) • 

1 L ( 50 * 20 > ♦ LP (20) 

COMMON /T I TL / Nl. TITLE l( l 5 ) • N 2 ♦ T I T LE 2 < 1 5 ) . N3 * T I TLE 3 ( 15) * NX, 

1 XTITLEI IS) .NHT.HORI Z(9) .NVT* VERT(9) .XHZ , Y HZ * THETA 1 * 

2 XVTi YV T. THET A2 .DCHAR ( 2 *2 > * NCHAR, CHFMT * EL A EEL • LLA BEL * 

3 DELTA, DELTA! .SUBDI V 
OATA CHOICE/*CRVFIT ■ / • SU BD I V / 1 » / 

DATA DELL* DEL 8 «PO, NL *RTYP. AQSC A » XO »XF* DE L T X , XF MT .RLMAX ,6MAX*RLMIN 

* /o* • 0*. 4* « C* 0* 1* #0**Q*.I* #2*0 *0* *0»* 0*/ 

DATA BM IN. YY.KOLOR /0 * * 1 0 C* , 1 / * EL A0 EL , L L AB EL/ 2 * 0/ 

DATA N1.N2*N3.XHZ.YHZ,THETA1.XVT,YVT,THETA2 /3*i,6*0./ 

DATA NX* XT I TLE/ 1 » 15* 1H / 

DATA NHT#HORIZ/l *9*1 h / 

OATA NVT.VERT /1.9*1H / 

OATA DCHAR/4*IH / . NC H A R/ 1 / . CHF MT /4 . 2 / 

DATA TITLE1*TI TLE 2 * T I T LE 3 /15*IH *15*1H «15*!H / 

OATA LP * BP /40*999*/ 

DATA DEL TA/IOO. / *D EL TA 1/10 0*/ 

END 

SUBROUTINE READ 

DI MENS ION A ( 4 > • C < 4)* IER( 3 > *U PLOT (20) *VFLOT(20) 

DIMENSION LABEL1 <9 ) • LA8EL2 < 9 ) 

REAL L • L P 

COMMON /VAL/ 0P ( 2 0 ) * B( 50 * 2 0 ) * F I ( 50 • 2 0 ) * F 2 ( 5 0 « 2 0 > • I GP < 2 0 ) * T ( 20 > * 

* L ( 50 .20) *LP ( 20 ) 

INTEGER RTYP* EL ABEL* RTYPEi UL ABEL • VLABEL 
COMMON / IN IT/ CELL .DELS *PO ■ NL , NF P * N I P * IDIM*RTYP 

COMMON /5CAL / CHO IC E * AB SC A * XO *X F , CELT X * XF M T * R LM A X * 0 MA X . RLM 1 N * 

1 BMIN.YY* KOLOR *NLP*NQP 
REAL *8 CHOICE 

COMMON / T I TL / Nl . T ITLE U 1 5 ) * N2 , T I TLE 2 ( 1 5 ) . N3 ,T I TLE 3< 1 5 ) , NX , 

X TI TLEI 15 ) ,NHT.HUR[ Z< 9 ) ,NVT.VERT(9) »XHZ *YHZ • THETA 1 , 
XVT, YVT* THET A2 ,DCHAR{2,2) . NCH AR , CHFMT * BL AH EL • LLA BEL* 
DELTA, DEL T A 1 ,SUBOI V 
DATA I ER/4 HRTYP*4H BL ,4HABSC/ 

NAMELIST /HEAD/ RT YP E. UMAX ,UMt N* VMAX • VMI N .KOLOfi, I P R * CHO I CE * PO « 

OELT AU , DEL T A V, DCHAR « NCHAR * CHFMT. UPLQT * NUPLCT , 
VPLOT* NVPLOT » A 65 C A , XO , XF *DEL TX , S UBC I V , XF MT , T ITLE I < 
N1,T ITLE2.N2 ,T I TLE 3, N 3 , X7I TLE, NX *L AB EL 1 • NL A B X , 

XPOS l* YPOS 1 * THET A 1 * L AO EL 2 , NL AB 2 , XPCS2, YPCS2 * 

THETA2 * ULABEL, VLABEL • X T ICK * YT [ CK 
EQUI VALENCE (RTYPE,RTYP) * ( D EL T AU . DELB ) ♦ ( CEL T A V * O ELL ) , 

( VFLOT (1 ) . LP< 1 ) ) . (UPLOT( 1 ) * BP < 1 ) ) • ( V « A X . RL M A X ) * 

( UMAX, 8MAX ) * < VM I N* RLMI N) * ( DM IN ,BMI N) • C NL AB 1 ,NHT ) * 
l LABEL 1 ( l > * HOR IZ ( 1))«(NLAB2. NV T) . ( LABE L 2 I 1 ) * VERT (1 ) ) *(XPGS1# XHZ) * 

< YPOS 1 , YHZ ) * (XP0S2 ,XVT) , (YFCS2iYVT) . (NVPLOT.NL P) . 

( NUPLOT » NBP ) •( UL ABEL *BL ABEL) , (VLABEL *L LABEL), 

7 ( XT ICK .DELTA). < YT I CK * OEL TA 1 ) 

READ ( 5 * HEAD ) 

IF ( RTYP * EG. 1) RTYP= 2 
LI = NLP + 1 
L2 ~ N BP F 1 
I F ( L l • GT • 20 ) GO TO 91 
DO 1 l =L 1.20 

1 LP( I > = 999. 

91 CONTINUE 

IF(L2.GT ,20 )GO TO 92 
DO 2 I=L2,20 

2 BP ( 1 ) - 999. 

92 CONTINUE 

IF (RTYP.LT.O.OR.RTYP* GT.2 ) CO TO 50 

IF ( A9SCA.LT *0 ••OR .AGSCA.GT. 1 * > GO TO 52 

IF < A 5 SC A* EG * 0* .AND* (XO* EG «XF) ) GO TO 62 

IF ( Nl *EQ.0.GR.N2*EQ*0 *OR. N3.EQ.0) WRITE<6. 700) 
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IF ( ( X HZ .L T • -1-25) .OR* (XHZ .GT. 8.C625)) WRITE(6.700> 

IF (<XVT *LT. -1.25) * OR • (XVT .GT. 8.C625)) WR1TE(6.700> 

IF ( ( Y HZ .LT. -.4) « OR • (YHZ • GT . 5.901) «J R I T E ( 6 * 7 0 0 ) 

IF < (YVT *LT. -.4) .OR* ( Y VT .GT. 5.90)) WRITE(6.7CQ> 

IF IKOLOR.LT. 1* OR .KOL OR *G T . 2 ) KGLOR = 1 

IF ( fcILABEL *LT #0 «0R .BLAOEL.GT *2 ) 0LA8EL = 0 
IF < LLABEL.LT .0 .OR .LLAEEL* GT .2 ) LLABEL - 0 
1=0 
NL = 0 
ML = 0 

IF (RTYP.EQ.l) GO TO 31 
READ FIRST FAMILY 
5 1 = 1 + 1 

READ < 5* 10 0) K ARDS • T ( I > 

IF (KARDS. EG.O) GO TO 30 
NL = NL + I 
DO 1C N=l* KAROS 
IFIRST = 4*(N-i)+L 

LAST = IFIRST ♦ 3 

READ (5.200) < B ( K * I ) ,F It K . I ) . K = I F I RST , L A ST ) 

10 CONTINUE 

DO 15 M= IF IRST *LAST 

IF ( 8<M* I> .NE.O..AND.F 1 (M. I) *NE*0. ) GO TO 15 
B(M. I ) = 999. 

FI (M* I ) = 999. 

15 CONTINUE 
GO TU 5 

30 IF (RTYP.EQ.Q) GO TO 99 

READ SECOND FAMILY IF NEEDED 

3 1 1 = 0 

32 1=1+1 

READ (5.100) KARDS.T(t) 

IF ( KARDS. EQ.O ) GO TO 99 
ML = ML + 1 

DO 35 N=l, KARDS 
IFIRST = 4* ( N— T ) + 1 
LAST = IF IRST 43 

35 READ (5. 200) ( L (K . 1 ) »F 2 { K » I ) .K = I FIRST *L AST ) 

DO 40 M— IFIRST .LAST 

IF ( L( M. 1 ) *NE • 0 * . AND *F 2 { M. l).NE«0.) GO TC 40 
L(M*1) — 999. 

F2 (M, I ) = 999. 

40 CONTINUE 
GO TO 32 
50 H= X 

GO TO 60 
52 N = 3 

60 WR1TE( 6.400) IER(N) 

GO TO 999 
62 WRITE( 6.600 ) 

999 STOP 
99 CONTINUE 

IF (ML.GT.NL) NL— ML 
IF ( I .NE .21 ) T< I > = 9 99. 

IFCRTYP .EQ.2 ) PET URN 
IF (IPR.EQ.l) CALL WRIT ( I PR ) 

WRITE( 6.801) 

IF ( I PR .NE * 1 ) WRITE (6.9 00 ) 

WRITEC6.800) T 
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100 
100 
200 
200 
20 1 
300 
40 0 


RETURN 

FORMAT !67X* I2*EX#F6.C) 

F ORM AT ( 67X * 12#€X«F5#2) 

FORMAT! 4{F8.S* IX* E 10 .3 ) • AX ) 

FORM AT ! 4! F6#5* E8.3# l X ) > 

FORMAT (2(F6#5»E8.3#16X),20X) 
FORMAT (12) 

FORMAT ( 1 HI . 31H********** VALUE 
TERM INATES. 


OF 


VARIABLE ,A4*52H NET IN CORREC 


IT RANGE. RUN 
268H4********* 

3 

480H********** 

5 

68flH********** 

7 

888H***A****** 

9 

500 FORMAT! 1 H 1 , 88 H ********* * 


// IX 


CORRECT RANGE 


0^2 


********** 

VAR I A 8LE 
******* ***/ | X . 

RTYP 

**********/ IX. 

NL 

**********/ 1 X , 

ABSCA 0 UR 1 

**********/ IX » I OH ** ********68X» l OH********** ) 

VARIABLES DELL AND DELB MIST BOTH BE 


2 — 2 0 


NO 


DESIRED. ALL FOUR VAR 


IN-ZERO. RUN TERMINATES. ********** 

600 FORM AT ! 1 HI # 99 H ******** ** IF LINEAR -SCALING 

| I A 3L ES XO.XF.DELTX.XFMT MUST Be READ. ******* / l X , 93H* ******* ♦ 

2 * CHECK TO MAKE SURE XO IS NCT EQUAL TC XF . 

3 ******* ) 

700 FORMAT ( 1HI . 1 0 2H ******** ** EARNING - GRAPH TITLES HAVE NOT BEEN S 

l PEC I FI EO. BLANKS U 1 LL BE SUBSTITUTED. ********** ///) 

900 FORMAT { 1 HI ) 

800 FORMAT! iHO. 10< 1X.F10 *3 » IX, )./. lX.10( 1X#F|0*3,1X) ) 

80 1 FORMAT ( IHO • 3X # 36HTHE-FQLL0 Wl NG-LINES-H AVE-BfcEN-I NPUT . > 

END 

SUBROUTINE WRIT (1FR) 

DIMENSION A! 4 > ,C< 4 ) , IER( 3 > . UPL QT (20) , VFLCT (20) 

DIMENSION L ABEL 1 (9 ) ♦ LA6EL2 ( 9 > 

COMMON /VAL/ BP(20 ) ,B( 50,20 *F 1 ( 50 .20) ( F2{ 50,2 0 ,IGP(20). T(20>. 

* L(50,20) #LP ( 2 0 ) 

INTEGER RTYP* EL ABEL# RT YPE* UL ABEL , V L ABEL 

COMMON / IN I T/ DELL .D EL B • PO .NL . NFP. N I P . I D I M , RTYP 

COMMON /SCAL / CHO 1C E. AE3SC A. XO .XF . DELT X . XFM T • RLM AX . BMA X . RLM l N. 

IBM IN . Y Y # KQLOR ♦ NLP# NBP 
RE AL *8 CHOICE 

COMMON /TITL / Nl,TITLEl(15).N2,TITLE2{ 1 5 > . N3 ♦ T I T L£ 3 l 1 5 ) . N X . 

1 XTITLE( IS > *NHT # HO RI Z( 9 ) »NVT,VERT(5) *XHZ * YHZ *THET A 1 . 

2 XVT, YVT.THETA2 .0 CH AR ( 2 # 2 ) * NCH AR . CHF M T . EL AB EL .LLABEL, 

3 DELTA# DELTA1 .SUBDIV 

EQU 1 VALENCE < RT YP E . R TY P ) . < DELT AU *DEL B > # ( CEL T AV .DELL ) . 

2 (VPLOT(l). LP { 1 ) ) *( UPLO T( 1) »6Pl 1) ). ( V N A X , RLM AX ) . 

3 ( UMAX. EM AX ) # ( VMI N* RLM I N ) • ( UM IN. BMI N) • ( NL AB 1 #NHT ) * 

4 ( L A B EL 1 ( 1 ) ■ HOP IZ ( 1)) *( NLA0 2.NVT) *(LABEl2(1) ,VERT(1) ) •< XPOS 1 . XHZ ) » 

5 ( YPOS1 * YHZ ) ♦ CXPO 52 #XVT ) • ( YFCS2# YVT ) , (NVPLOT .NLP ) # 

e> ( NUPLOT . NBP ) . < UL A 8 EL .BLA0EL ) .{VLAOEL .LLABEL) . 

7 { X T ICK . D EL TA ) . <Y TICK, DELT A 1 ) 

WR I T E ( 6 • 1 ) IPR 
WRIT£(6,21 
WRITE(6,3>RTYPE 
WRITE! 6.4) PO 
WRITE( 6. 5) 

WRITE! 6.6) 

WRITE! 6#7) 

WRITE! 6. 8) 

WRIT E( 6*9) 

WRITE! 6. 10 ) 


DEL T AU 
DELTAV 

(VPLCT! l> .1=1 .NVPLCT) 
NVPLOT 

( UPLGT ( I > #1=1 .NUPLOT) 
NUPLOT 


WRITE! 6. 11 ) 

MR ITE! 6. 12 > 
WRITE! 6*13) 


CHO ICE 


SI 



WRITE<6.14> UM/IX 
WR ITE* 6. 15) UM IN 
WRITE*6.16) V M A X 
WRITE<6.17) VMIN 

wRire(6.ia> kolor 

WRITE* 6* 19) ABSC A 
WRITE* 6* 20 > XO 
WRITE(6.2l> XF 
WR I T £{ 6 * 22 I DEL T X 
WR ITE* 6* 23) SL80 I V 
WRITE(6*24> 

WRITE! 6* 25) < T ITLE l* t > . 1 = 1 • 1 S) 

WRITE* 6 » 26 ) N 1 

WRITE* 6. 27) (TITLES* II *1=1 .IS) 

WRITE* 6*28) N 2 

WRITE* 6,29 ) ( T ITLE 3* I) * 1 = 1 .15) 

WRITE*6,30> N 3 

WRITE* 6* 31 ) (XT ITLE* 11.1 = 1*151 
WRITE(6.52) NX 

WRITE* 6. 32 > < L ABEL 1 ( I) .1 = 1 .9) 

WR I T E ( 6 . 33 1 NLAQ1 

WRITE<6.34> XFCSl 

WR ITE* 6. 35 ) YPCS1 

WRITE{6. 36 ) THETA I 

WRITE* 6. 37) C LABELS* 11.1 = 1.9) 

WR I T E ( 6 « 38 ) NL AS 2 
WRITE* 6. 39 ) XPCS2 
WRITE* 6. 40 ) YP052 

WRITE* 6.41 ) THE T A2 
WRITE* 6.42 ) ULABEL 
WR ITE*6.43 ) VLABEL 
WK ITE* 6*44 ) XTICK 
WRITE(6.45> YTICK 

WRITE* 6. 46) C ( OCHA R (, 1 . J ) .1=1. 2). J= 1.2) 

WRITE* 6.47 ) NCFAR 
WR ITE* 6 • 48 ) CHFMT 
WRITE* 6. 49 ) XFMT 
RETURN 

1 FORM AT ( 1HI .5X .6FIPR ® .12) 

2 FORMAT* 1H0.3QX.33HCARPET PLOT CONTROL PARAMETERS) 

3 FORMAT ( 1H0 . 5X * EHRTYPE = .11) 

4 FQ RM AT ( I H .5X.SHP0 = .F4.1) 

5 FORMAT { 1 H .5X.SHDELATU = .F8.5) 

6 FORMAT* 1H # 5 X . 9HDELT AV = .F8.5) 

7 F ORM AT * 1 H . 5 X . EHV PLO T = , / .5 ( 1 PE 1 0.3 . 3 X ) ) 

8 FORMAT* 1H .5X.9HNVPL0T = .12) 

9 FORMAT* IH .SX.EHUPLnT = ./ . 5 « l PE 1 0 . 3 . 3X > ) 

10 FORMAT* 1H , 5X , 9HNUFL CT = .12) 

1 1 FORMAT* IHO. 30 X.25H INTERPOL AT ION PARAMETERS) 

12 FORMAT* 1H0.5X.9HCHCI CE = . A8 ) 

13 FORM AT* IHO. 30X .16HGR ID PARAMETERS) 

14 FORMAT* IHO. 5X.7FUMAX = ,F4.1) 

15 FORM AT ( l H .5X.7FUMIN = *F4.1> 

16 form at < i h .sx.thvmax = .F4.i) 

17 FORMAT ( 1 H .SX.7HVMIN = .F4.1) 

18 FORMAT* IH .5X.EHKOLOR = .It) 

19 FORM AT < IH .5X.6FABSCA = .F4.1) 

20 FORM AT * 1 H .5X.SHX0 = , 1PE12.E) 

21 FORMAT* IH .5X.5HXF = . IPE12.E) 

22 FORMAT * 1 H .5X.CHDELTX = .IPE12.5) 

23 FORMAT* IH .5X.9HSU90IV * .1PE12.5) 

24 FORMAT* IHO. 30 X , 34HLA BEL I NG AND TITLING PARAMETERS) 

25 FORMAT* IHO. 5X.SHTITLE1 = .15A4) 

26 FORMAT* IH .5X.SHNI = .12) 
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27 

FORMAT ! 1H 

. 5X.SHT ITLE2 

- . I 5 A 4 ) 

28 

FORMAT ( 1H 

.5X.5FN2 = • 

12) 

29 

FORMAT ( 1H 

• 5X ♦ 9HT ITLE3 

= . I5A4 ) 

30 

FORMAT! 1 H 

.5X.SHN3 = . 

12) 

31 

FORMAT ( 1 H 

.SX.9HXTI TLE 

- . 15A4 ) 

52 

FORMAT { 1H 

.5X.5HNX = . 

12) 

32 

FORMAT! 1H 

. 5X » 5HL ABEL 1 

- « 9A4 ) 

33 

FORMAT ( 1 H 

• 5X * 6HNL AQ 1 

= .12) 

34 

FORMAT t 1H 

« 5X . 6HXPCS 1 

= ,F8.E) 

35 

FORMAT! 1H 

.5X.EHYPCS1 

= *F 6 • E ) 

36 

FORMAT! 1H 

« 5X*9FiTHET A1 

- ,FS.2> 

37 

FORMAT ( 1 H 

. 5 X . 7HLA B2 = 

. 9 A 4 ) 

38 

FORMAT! 1H 

.5X.EHNLAB2 

= .12) 

39 

FORMAT! 1H 

»SX .EHXPCS2 

- »F8.5) 

40 

FORM AT < 1 H 

•5X| €HYPCS2 

= »F8« E) 

41 

FORMAT ( I H 

.5X.9HTHETA2 

= , F8 .5 ) 

42 

FORMAT ( 1H 

.5X .ShULABEL 

= .11) 

43 

FORMAT { 1H 

* 5X , 9FVL ABEL 

• . 11 ) 

44 

FORMAT! IH 

. 5X * 6HX T ICK 

- .F6.2) 

45 

FORMAT C 1H 

,5X . EH YT ICK 

= »F 6 • 2 ) 

46 

FORMAT ( 1 H 

♦ 5X. EHDCFAR 

= ,2A4) 

47 

FORMAT! 1H 

* 5X . EHNCHAR 

= .11) 

48 

FORMAT ! IH 

, 5X. EHCHFMT 

= ♦F8.E) 

49 

FORMAT! IH 

*5X » 7hXF MT = 

*F8 .5 ) 


END 

SUBROUTINE SUPER ! XMIN. XMAX tYMIN* YMAX ) 

COMMON/PLQ T/AC 200 ) 

DIMENSION PVX{ 2 ) .PVY(2 ) 

ADDS REGISTRATION MARKS 

XM IN.X MAX, YM I N AND YMAX ARE RASTER COCOS OF GRID* 

X 1 =X M I N— 35 0 » 

X2=XMAXf3SO. 

Y l=YMl N-100* 

Y2=YMAX+10O« 

CALL 06JCTG! A.X I* Yl, X2.Y2) 

CALL SETSMGI A » 14.1 • ) 

I T I ME=0 

CALL S ET SM G! A * 3 0 « 0 «5 ) 

10 I T IME= IT IME+l 

GO TO ( 1 .2 * 3*4)* I TIME 

1 XC=250* 

YC=YMIN 
GO TO 20 

2 XC-250. 

YC = Y MA X 
GO TO 20 

3 XC-3850* 

YC=YMAX 
GO TO 20 

4 XC-3850. 

YC-YMIN 

20 PVX< l)=XC-75.0 
PVX ( 2)~XCF75« € 

PVY(1)=YC 
PV Y { 2 > =Y C 

CALL L I NESG( A . 2»PVX»PVV ) 

PVX( 1 > =xc 

PVX( 2) =xc 

PVY< l)=YC“75.0 

PV Y ( 2 ) =YC + 75 • 0 

CALL L INESGCA. 2 iPVXt PV Y ) 

CALL 5£TSMG< A • 03. 5*0 > 


53 



CALL C IR ARG( A* AC *YC*37.5*0«*36C* ) 

I F ( l T I MB .LT • 4 ) GO TO 10 
CALL SET SM G ( A * 30 « 1 • ) 

RETURN 

END 

SUBROUTINE IN I T AL 
INITIALIZE NAMELIST 
REAL L*LP 

COMMON / V AL/" EP ( 20 ) • 8 ( 50 * 2 0 ) »F 1 1 5 0 * 2 0 ) *F 2 I 5 0 *2 C ) * I GP (2 0 > * T ( 20 ) 
1 L(50*20), LP (20) 

COMMON /PSI/ PSI 1(50*20) ,PS12(5C*20> 

COMMON /LABS/ L A8 1 ( 8 ) * L A 62 ( 8 ) 

DATA I BLANK/ 1 H ✓ 

X D 1 M = 50 
DO 50 1=1*8 

LABI ( 1 j = IBLANK 
LAB2(I) = 16LANK 
SO CONTINUE 

OO 100 1=1*20 

IGP< I > = 0 
DO 100 JslilDIM 
PSi 1 ( J * I > = 999* 

PSI2< J * 1 ) = 999. 

B( J* I ) = 999. 

FI ( J. I ) = 999 . 

F2(J.I) = 999* 

100 L( J, I ) = 999. 

RETURN 

END 

SUBROUTINE GRICPT ( RP * A * I GP ) 

DIME NS ION RP( 20) *A (50, 20 > . 1GP< 20 I 
INTEGER RTYP 

COMMON / I N I T/ CELL • DELB ,POiNL » NF P, N I P * IOIM*RTYP 
DO 40 1=1* NL 
P = RP ( I ) 

IF tP.NE.999.) GO TO 20 
GO TO 99 

20 OO 30 J=1 , 10 1M 

IF ( AC J* 1 ) *NE«P) GO TO 30 
I GP ( I ) = J 
GO TO 40 
30 CONTINUE 
40 CONTINUE 
99 RETURN 
END 

SUBRQUT INE SCALE (PSI * RNUH.PHI) 

DIMENSION PS I; (50* 20) 

REAL L.LP 

COMMON /VAL/ ep<20 ),B< 50,20) ,F 1( 50,20) ,F2(S0,2C> *IGP(20) ,T<20> 
1 L(50* 20 )• LP( 20) 

INTEGER RTYP 

COMMON / I N I T/ CELL, OELB.PO *NL*NFP, NIP , 1DIM.RTYP 

SCALES CARPET FLOT IN VERT DIRECTION. 

DO 110 1 = 1 * N I P 

P = PS I ( 1* I ) 

IF ( P . EG * 999 • ) GO TO 110 

PMN = P 

PM* = PMN 

IMN = I 

KMN = 1 
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GO TO I 1 1 

110 CONTINUE 

111 DO 115 1 = 1 »N IP 
00 115 J=1 i NFP 
K a IGP(J) 

PP = PSIiK, I ) 

IF (PP.EQ.999.1 GO TO 115 

IF (PP.GE.PMN) GO TO 114 

PMN -PP 

KMN - J 

1MN = 1 

GO TO 115 

114 IF i PR • GT • PMX 1 PM X=PP 
US CONTINUE 

RNUM = I NT ( 2 * *PQFP MX —PMN ♦•01 J 
IF {RTYM.eo.O) GO TO 120 
ITEMP = IMN 
IMN = KMN 
KMN - ITEMP 
120 AMI NL = LP(KMN) 

BMIN = BP (IMN) 

PHI = PO - AMINL/DELL - BMIN/DEL6 

RETURN 

ENO 

SUSROUT I NE LlKT (AiTtF) 

DIMENSION A( 50 ,20 ) ,T<20) ,F f SCI, 20 > 

INTEGER RTYP 

COMMON / I N I T/ DELL ,DELB ,PQ ,NL ,NFP,NI P, ID IM, R7YP 

INT a IOIM/NL 

IL = NL — 1 

DO 20 1=1, IL 

K 1 = I NT 4 { I— 1 ) F 1 

K2 = Kl F INT - 1 

K = K1 + l 

KN - K2 F 1 

A ( K I « L ) = TCI) 

DELA = T{ I Fl ) - TO) 

DEL— DEL A/ I NT 
DO 20 I T= 1 » N IP 
FI = F( K 1 • IT) 

F2 - F(KN, IT > 

IF C ( F 1 .EU.999 « ,OR« F2 *EQ • 999. ) • AND • I T *NE . 1 > GO TO 20 
FAC =* { F 2— F l ) / CEL A 
DO 10 N=K * K 2 

IF (IT.EQ.l) AIN, 1 ) = A ( N— 1,1) F DEL 
IF C FI .EQ.999 • .OR .F2 ,EQ ,999. ) GO TO 10 
F t N , IT J ^ FI F ( A ( N » 1 ) - A (K 1# 1) ) *FAC 

10 CONTINUE 
20 CONTINUE 

1ND = INT4IL F 1 
A ( IND, 1 > = T(NL 1 
RETURN 
END f 

SUBROUTINE REVERS I X, Y,LGP) 

01 MENS ION X( 50,20) *Y(50,20), XT(50) ,YTC £0) ,LGP< 20 ) 

REAL L »LP 

COMMON /VAL/ EP ( 20 > , B ( S 0 * 2 C ) ,F l( 50 . 20 ) ,F 2 ( 5 C ,2 C > , l GP (2 0 > , T ( 2 0 > , 
1 L < 50, 20 ) • LP< 2 0) 

COMMON / I N 1 T / DELL, 0ELB,PC, NL ,NFP, NIP, IDIM.RTYP 
DO 25 1=1. NFP 
J = NFP - I F 1 
25 LGP(l) = IDIM - I GP ( J ) F l 
DO 100 1=1, NIP 
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DO 50 J~1 t ID IM 

K=IDIM-J*1 

XT ( J ) = X ( K * 1 ) 

50 YT < J J = Y(K f I ) 

DO 100 N=1 * I D IM 
X< N* I ) = XT! N> 

100 Y ( N« I ) - YT(N > 

RETURN 

END 

SUBROUTINE FINDF!FP*A* AY»CY* T) 

REAL L * LP* LL 
INTEGER RTYP 

COMMON /IN IT/ CELL .DELB. PO.NL.NFP, NIP, IDIM, RTYP 
DIMENSION FP120)tA(50# 20)*AY!5G# 20) *CY(5C*2Q) * T ! 20 ) 
DIMENSION LAB EL (2) •MIS(20) 

DATA LABEL/4H X2=*4H Xl=/ 

IF (RTYP.EQ*1) GO TO 5 

I TP = LABEL ( 2 ) 

IRD - LABEL < 1 ) 

LTYP = 2 

GO TO 6 

5 I TP = LABEL! I ) 

I RO = LABEL! 2 > 

LTYP= 1 

6 DO 7 I^liNIP 

7 mis( n = o 

LIM = IDIM * 1 
INT - IDIM/NL 
DO 50 IL= 1 , NL 
IP - INT*< IL-1 ) + l 
DO 20 1-UNIP 
BFIX s FP< 1 ) 

DO 15 J^l.LIM 
PI = AOtlL) 

IF { PI .£Q#9<39# > GO TO 14 

1N2 - J+l 

P2 “ A { IN2, IL ) 

IF < P2.EQ.999. ) GO TO 14 
XI * AY( J* IL) 

IF (PI# NE . 8F IX) GO TO 10 
CY ( IP* I > ~ xi 
GO TO 20 

9 CY! IP. I I - AY( IN2 • IL ) 

GO TO 20 
10 CONTINUE 

IF ( P2 • LT • BF 1 X J GO TO 15 

IF { PI .LT.BFIX.AND.P2.GT #BFI X ) GO TO 16 
IF ( P2 • EG* BF IX) GO TO -3 

15 CONTINUE 
14 CONTINUE 

MIS! I J = MIS! I ) + 1 
IF (P2.LT*6FIX) GO TO 50 
GO TO 20 

16 CONTINUE 

X2 = AY( IN2* IL ) 

CY ! I P» I ) = XI + ( BF 1 X— P 1 I * { X2— X 1 1/(P2— FI ) 

20 CONTINUE 
50 CONTINUE 
IER = 0 

II =1 

55 DO 60 1=11, NIP 
IN = I ♦ IER 

IF (MIS! IN ) #LT#NL1 GO TO 60 
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WR ITE(6«100l> ITP.FPU >«LTYP 
IER = IER + 1 

J =1 + 1 
GO TO 70 
60 CONTINUE 
GO TO 999 
70 DO 80 N^J.NIP 
80 FP(N-1 ) — F P ( N ) 

FP(NIP>= 999* » 

NIP = NIP - 1 
II - J-l 

IF (II *LE* NIP) GO TO 55 

WRITE! 6. 1002) NIP. (FP( II). 11=1. NIP) 

999 RE TURN 

1001 FORM AT { l Hi * 9H THE V AL UE . A 4 • F6 .2 . 4 OH DOES NOT FALL OK ANY LINE UF CO 
INSTANT X.IU1H.) 

1002 FORMAT !5X. 43HTFE NUMBER OF I NT EH POL A TEC LINES CHANGED T0.I3.30H.TH 
IE NEW ARRAY IS AS FOLLOWS -/5X.2GF6.2) 

END 

SUDROUT INE INTERP < Xl N. Y IN. LGF . INT.M.N. IERR ) 

SPLINE OR CRVFIT I NT ER FLLA T I ON DETERMINED BY CHOICE 

COMMON /SCAL / CHO ICE. A BSC A. XQ.XF.DELT X. XFMT . RLM A X .BMAX.RLM IN. 

1 BM IN .Y Y » KOLQR .NLP. NBP 
RE AL *8 CHOICE 
REAL *8 SPLINE 
DATA SPL INE/8FSPL INE / 

DIMENSION XiN(M.N) .Y IN(M.N ) * LGP( 1 ) 

COMMON X { 50 ) * U ( 50 ) ,S(5C)*0EL(50) .Ql.ON.NFTS 
I OE0UG=O 

IF (I£RR *EQ. C) GC TO 150 

I D£BUG= 1 

WR1 T £ { 6* 100 ) 

100 FORMAT ( 1H1 . / / .4 9X . 35H + + 4 44 INTERP HAS BEEN ENTERED *44*4,//) 

WRITE! 6* 110) M.N.INT 

110 FORM AT(9X*6H H — .I6.9X.6H N = . 1 6 * 1 1 X . 6H [ NT — .16.///) 

WRITE! 6. 1 15) <LGP< I > .1=1 ,N ) 

115 FORMAT (A- OX . 

154HTHE INDEX OF THE FIRST PLOTTED POINT FOR EACH CURVE IS.//. 

24! 40X.5! I2#8X). />•///> 

DO 140 J = 1*N 
WRITE< 6. 120 ) J 

120 FORM AT ( l X. // * 4 EX. 1 2HCUR VE NO. = , I 6 , 5X . 1 3HX I N. Y l N FAIRS*//) 

WRITE(6. 130) (XIN( I • J ) .YlN(I.J). 1 = 1. M) 

130 F ORMATI1X.41F10.2.2X.E10.2.6X) ./) 

140 CONTINUE 
150 CONTINUE 
NCUR VE=0 
I ERR=0 
01 = 0 .0 
QN=0.0 
Q 1 STRT-=0 « 0 
QN ST RT = 0 . 0 

IF ( ( INT .GT. 0) .AND. (INT .LE. M)) GC TO 1500 
I ERH=1 
GO TO 9999 
1500 CONTINUE 
2000 NCURVE=NCURVE+ I 

IF INCURVE • GT « N) GC TO 9999 
DO 2100 XFIRST=1*M 

IF IXIN! IF IRST.NCU RVE) .NE. 999.) GO TO 30C0 
2100 CONTINUE 

GO TO 9999 
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3000 NPT5— 1 

X< 1 >=X IN( IF IRST «NCURVE) 
U( 1 )=YIN( 1FIPSUNCURVE) 
OO 3200 I=IFIRST,M 




IF ( X I N C I • NCUR VE ) .EG. 999.) GO TO 3600 



IL AST= I 


3200 


CONTINUE 


3600 

DU 

3620 1= UN 




IF (LGPCI) .LE. 0) GO TO 

3 G 5 0 



IF CLGPC I ) • GE . IF IRST ) 

GCJ TO 3640 

3620 


CONT INUE 


3640 

IF 

CLGPC I) ♦LE. ILAST) GO TO 

3700 

3650 

IERR=2 



GO 

TO 9999 


3700 

INDEX=LGP< I ) 



IF 

( INDEX .EG. IF IRST) GO TO 

4000 

3000 

IF 

I { I NDEX — INT) .LE. IFIRST) 

GO TO 4000 


INDEX= INDEX- INT 



GO 

TO 3800 


4000 

IF 

(INDEX • €E • ILAST) GO TC 

4200 


IF 

(INDEX .EQ. IF IRS T ) NPTS= 

NPTS-l 


NPT S=NPT SF 1 

XCNPTS )“X|N{ INDEX, NCUR VE > 
U( NPTS )=YINl INDEX, NC UR VE ) 
INOEX= INDEX* INT 



GO 

TO 4000 


4200 

NPTS=NPTS+ 1 

XCNPTS ) =X I N{ ILAST, NCUR VE } 



U ( NPT S ) = Y I NC I LAST. NCUR VE > 

IF {CHOICE .NE. SPLINE) 60 TO 6000 
IF ( A6SC A ♦ EQ ♦ 0 • ) GO TO 4400 
DO 4300 I=1*NPTS 

U( I ) -AL06 1 0 { U { I ) ) 

4300 CONTINUE 

4400 IF ( NP T S .L E « 3 ) GO TO 5400 
Ql-{ U( 2)-U< 1 ) )/( X< 2) -X< l > > 

QN={ U{ NPTS)~U< NPTS-1 ) ) /( X< NPTS >-X< NPTS- I ) ) 

G 1 STRT=QI 
GNSTRT=QN 
CALL SPLN2 
NUMPT = I 

5000 DO 5200 I— IFIRST* ILAST 

IF (XINC I «NCURVE ) .NE. X<NUMPT>) GC TO 5100 
NUMPT = NUMP T ♦ l 
GO TO 5200 

510 0 Y I N ( I • NC UR VE ) = SPLI N< XI N( UNCURVE) ) 

IF ( ABSCA * E Q ♦ 1.) Y I N ( I , NCURVE ) = I 0 . 4 * ( Y I N ( I . NCU RVE ) ) 

5200 CONTINUE 

5300 CONTINUE 

GO TO 5800 

5400 IF ( NPTS «EQ« 21 GQ TO 5700 

A = ( { X( 1 >-XC 2 2)-UC3 ) > - ( X{ 2>-X< 3) ) *<UC 1 > -U<2) > ) / 

I i (X( 2)**2-X ( 3 )**2) *C X( 1 >-X <2>) - ( X{ 1 1 ** 2-X ( 2 > **2 ) *C X { 2 ) -X ( 3 ) ) ) 

6 a { U(l) - U ( 2 ) - A*( X(l)**2. - X(2)**2. ) ) /< X(l) - XC2) ) 

C a U( 3 ) - A*X{3)**2. - B+X{3> 

MLAST = ILAST - INT 
DO 5600 1= IF [RST, MLA$T * INT 
11= I + l 
12=1 + INT - I 
00 5600 ND=I 1 , 12 

Y I N { NO » NCU RV E 1 = A *X IN ( ND • NCUR VE ) *# 2 • ♦ 8 * X I N ( NO * NC URV E ) + C 

5600 IF ( ABSCA .E Q. 1 ♦ ) Y I NC ND . NCUR VE I = 1 0 . ** Y I N ( ND » NCU H VE ) 

GO TO 5800 
5700 CONTINUE 
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ND1 = IFIRST + 1 

ND2 - IL AST - 1 

DO 5750 I = N0 1 < ND2 

Y l N! I* NCURVE ) = U{ 1) + C X l N ( I • NC UR VE ) - X ( 1) > / C X ! 2 ) - X ( 1 > ) 4 
1 (U(2) — U < 1 ) ) 

5750 IF ( ABSCA « EQ « 1 • ) Y IN ( I « NCU RV E > = l 0* 44 Y I N < I * NCUR VE ) 

5800 IF C 1DE0UG *EG* 0) GO TO 2000 
WRITE! 6*200) 

200 FORMAT ( IHI * // *48X* 27H44444 INTERMEDIATE INTERP DEBUG 44444,//) 
WfUTE!6*210) NCURVE* NPTS* I F l RST* ILAST 
210 FORMAT ( 9X, 8HNCURVE = * I 6 * 5X *6HN PTS = • I 6 » 5 X # 8H I F I R S T =,I6*5X* 

I 7 H IL AST = • I 8 * // ) 

IF ( NP T S T * 2 ) WRITE<6*220 > 

220 FORM AT ( 1 X* /// * 

1 45X ♦ 4 4H444*44444*4**4444444444 44 4444444444444444444*/, 

2 45X.44H4 THIS CURVE HAS FEWER THAN TWO POINTS 
345X*44H4 NO INTERPOLATION HAS BEEN CONE ♦*/. 

4 45 X* 44 H****** ******* ****** 44 4 44444444444444*4444444 * ///) 

WRITE! 6, 23 0 I C 1STRT* QNSTRT *Q 1 » QN 
230 FORMAT!//, 57X, 11HFIRST PG I NT , I CX * 1 OHL A S T POINT*//, 

146X#6HSTART ,eX,F9.4*llX,F9.4,//* 

246X,6HFINISH,€X*F9.4,l 1X*F9.4,// ) 

WRITE! 6, 120 ) NCURVE 

WR I TE( 6, 13 01 ( X I N < I * NCURVE ) * VI M I • NCURVE ) * I - 1 « M ) 

GO TO 2000 
6000 CONT INUE 

IF < ABSCA • EO • 0.) GO TO 6200 
DO 6100 I=1,NPTS 

U( I ) = ALOG 1 0 ! U ( U > 

6100 CONTINUE 

6200 CONTINUE 
M0Dt=l 

ii=npts 

12 = 1 NT 

CALL CRVFIT! MODE, I I • X* U * 12 , l 3 * X I N< IF I R S T , NC UR V E J , Y I N C I F I RS T , NC UR VE 

*) 1 

IF (ABSCA * £ Q • 0#) GO TO 6E00 

DO 6400 |= IF IRST, ILAST 

YIN! I * NC URVE)— 10 • 4 4 ( Y 1 N< I*NCURVE>) 

6400 CONTINUE 

6500 CONTINUE 

GO TO 5800 

9999 IF ( IDE8UG • EG • 0) GO TO 7000 
WRITE! 6, 300) 

300 FORMAT! IHI *// *EOX, 34H4444* INTERP IS BEING EXITED 44444*//) 

WRITE! 6* 3X0) I ERR* NCURVE, I F I RST* ILAST, KPT S 
310 FORMAT ( 9X, 6H IEFR = , 1 6 * 5 X * 8HNCUR V E = * I 6 * 5X , 6H IF I R ST =*I6*5X, 

1 7HILAST = * 16, 5X, 6HNPTS =,I6,//) 

WRITE! 6. 110) M*N. I NT 
DO 340 J=1*N 
WRITE! 6* 120 ) J 

WRITE (6*130) ! X IN ( I ♦ J ) *YIN{1*J1*I=1*M) 

340 CONTINUE 

400 FORMAT ( SX* 2HA= ,F6* 2* 3H B=*F6*2,3H C=*F6,2) 

7000 RETURN 
END 

SUBROUTINE CR V F I T ( MOCE , L 0 , X,Y*MO*NO*U,V) 

INTERPOLATION EY CRVFIT 


DI MENS ION 
DIMENSION 
EQUIVALENCE 

1 

EQUIVALENCE 


XC50>* Y( SO > * U( 50 ). V(50) 

A0(2)*B0(2) 

(A, AO! 1) )#(0*AO(2) *00! 1) )»{C*BCC2> ) • 
(P0*X2). (QO* Y2) *(DX*A2) * (DY*B2) 
<FLM*TS,Z ) * C JP * JS) • 
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x (DU* DA *D »X1) ,(DV .0B*R*Y1 ) * 

2 <S2*S20*A1 ) i (S3,S03*B1 ) • 

3 (PL.S12) »CP2«CI2l*(P3»Rl2) * 

4 ( Ql» SI 2) *(Q2*C13)» (Q3,R13) 

SCR(SIJ.CI J > = AGS(S IJ )-ABS(CI J>*1 .OE-8 

MD=MODE 

L = LQ 

LPS=LO 

M=M0 

KP1=L*M+1 

1P-L+1 

DO 10 J P= I * L 
KPi =KP1 -M 
IP- I P -1 
U(KPU-X( IP > 

V t KP 1 > = Y< IP ) 

10 CONTINUE 

KP2- 1 
KP3= 1 

OO 20 1=2, L 

KP2-KP2+M 

IF< U ( KP 2)*EG,U( KP 31 #AN0*V(KP2 ) #EQ • V (KP3 J ) GO TC 20 

KP3=KP3*M 

U(KP3)=U(KP£) 

V<KP3)-V(KP2) 

20 CONTINUE 

L=KP3/M*1 
N=KP3 

IF(N.EQil) GO TO 690 
I F { MO * NE m2 ) GO TO SO 
30 OO 40 KP4=l*N,M 
TS=U(KP4) 

U(KP4)=VIKF4> 

V (KP4 )=TS 
40 CONTINUE 

50 MM1=M-1 

DZ=1 ,0/M 

iF<L.EQ,2) GO TO 100 

LM1=L-1 

GO TO 200 

100 DU=(U< N)-U{ 1 } )*DZ 
DV=t V( N 1 —V ( 1 1 )*0Z 
00 110 KS=1*MM1 

U(KS+ll=U<KS)+DU 
V(K$+1 )=V«KS)+DV 
110 CONTINUE 

GO TO 800 
200 X3=U<1) 

Y3=V( 1 ) 

X4=U(M+1 > 

Y4=V(M+1> 

K 5— l fM+M 
X5=UtK5l 
Y5=V(K5) 

A3=X4— X 3 

B3 = Y 4— Y 3 

A4=X5-X4 

04=Y 5 — Y4 

534=A3*E4- A4 + 83 

C34=A3*A4+B3*B4 

R34=SCR(S34,C34) 

1 F ( MO* LE , 3 ) GC TO 400 

K1=N-M-M 

X1=U(KI) 
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Yl^VCKI I 
K2=KlfM 
X2 — U (K2) 

Y2= V < K 2 ) 

A 1 — X 2— X 1 
82— Y 2 - Y 1 
A2=X3-X2 
62-Y3-Y2 

230 S12=At *B2-A2*B1 
Cl 2=A1 *A2*Bl *62 
Ri2=SCRI S12.C 12) 

S I 3= A 1 ♦8 3'- A3 *0 1 
C 13= A1 *A3*B1 *63 
R13=SCR(S13*C13) 
S23=A2*B3-A3*B2 
C2 3=A2 *A3+B2*B2 
R2 3= SCR ( S23. C23 > 
S24= A2*B4- A4*B2 
C24=A2* A4+B2*B4 
R24-SCR< S24* C24) 
ASSIGN 240 TO LBL 
GO TO 500 

240 DO 290 I = 2*L 
KS=K5+M 
X2=X3 
Y2 = Y3 
X 3=X 4 
Y3-Y4 
X4-X5 
Y4 = Y5 
A2 = A3 
B2 = B3 
A3 = A4 
B3 = B4 
S 1 2— S2 3 
R I 2=R23 
SI 3=S24 
R 1 3=R24 
S2 3- S3 4 
C2 3— C34 
R2 3— R3 4 
COS2=SGN*COS3 
SIN2=SGN*S1N3 


IFC ULT.LM1) 

GO 

TO 

270 

IF <MD*£Q.4 > 

GO 

TO 

26 0 

1FI I *EQ.LM 1 ) 

GO 

TO 

45 C 

A4-A5 
t34"=B5 
GO TO 260 
IFIK5.GT.N) 

K ^ = 

1 + M 



270 X5=U(K5) 

Y5=V(K5> 

A4=X5-X4 
84 - Y 5- Y4 

260 S24= A2 *B4- A4*B2 
C24-A2*A4*B2*64 
R24-SCRC S2 4 * C 2 4 ) 

S 34= A3 * 04— A 4 *63 
C34-A3*A4+B3*64 
R34=SCG( S34.C34) 
ASSIGN 600 TO LBL 
GO TO 500 
290 CONTINUE 
GO TO aoo 
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400 IF ( R34 «LE« 0* C ) GO TO 430 
IF(MO.EQ.3) GO TO 420 

410 A 1 — A 3 
A2=A3 

0B=2 «* A3*S34/< A 4*1 A3+A4 ) ) 
02=03-06 
61=62-00 
GO TO 230 
420 0A=A4~A3 
00 = 8 4— B3 
A2=A3-0A , 

62=83-06 
A 1 = A2— 0 A 
61= 02 — Oft 
GO TO 230 
430 A 1 = A 3 
A2 = A3 
01=03 
G2 = B3 
GO TO 230 

450 IF (R23.LE.0.0 1 GO TO 480 
I F < MD« EQ * 3 ) GO TO 470 

460 A5=A3 
A4= A 3 

L>0 = 2 • *A3*S23/ r ( A 2*1 A2+A2) ) 
84=03+00 
05=04+08 
GO TO 280 
470 DA=A3-A2 
DB=B3-B 2 
A4 = A 3+ O A 
B4=S3+DB 
A5=A4+DA 
05=84+ D0 
GO TO 280 
480 A5=A3 
A4 = A3 
65=03 
04 = 83 
GQ TO 280 
500 5GN=1.0 


IF1R23.LE.0.0 > 


GC 

TO 

550 

l F ( (R1 2.LE.0 .0) .AND. (R34.LE.0.C) 

) 

GO 

TO 

580 

IF ( < R13.LE.0.0) .AND. (R24.LE.0. 0) 

> 

GO 

TO 

580 

IF ( { R 1 2.LE.0.0 ) .OR • 1RZ4.LE.C.0) 

) 

GO 

TC 

560 

IF ( ( R13.LE.0.0) -OR- (R34.LE.0.0) 

) 

GO 

TO 

570 


52=S 1 2*524 
S3=S 13*534 
PAT CH= 5 2 48 3 

IF1PATCH.LT.0.0 ) S 3 =— S 3 
A= S2 * A 3 * A3 — S 3* A 2* A 2 
8= S2*A 3*03 —5 3* A 2* 8 2 
C=52*B3*B3-S3*E2*B2 
0= S2 3* SORT 1 S2 * S3 ) 

P A T C H= B * 0 

IFIPATCH.LT #0*0 ) D=-D 
8=B+D 

S 2 0= A2 * B 0 ( l >-AC< 1) 4B2 
S 0 3= AO 1 1 ) * 03— A3*B 0 ( 1 ) 
PATCH=S20*S03 

IFCPAT CH.L6.0.0) GO TO 510 
COS3=AO( 1 ) 

SIN3=801 1 > 

GO TO 520 
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blO S20- A2*B0{ 2 )-AO{ 2 ) *B2 
C0S3=A0<2) 

SIN3-80 (2i 
520 PATCH=S20*S23 

I F ( PAT CH»GT *0*0) GO TO 59C 

C0S3=-C0S3 

SIN3--SIN3 

GO TO 590 

550 IF ( C23 «LT * 0 * C ) SGN--!*0 

560 COS3=A2 
S|N3=B2 
GO TO 590 
570 C053-A3 
S I N3 — B3 
GO TO 590 
580 CQS3=A2+A3 
S IN3=B2+83 

590 lF(«0.Lfc*2) GO TO L BL * C 24 C ♦ 6 0 0 ) 
R— SQRT < COS34C0S3+S IN3*SIN3 ) 
COS3=COS3/R 
SI N3-SIN3/R 
GO TO LOL* (24C t 60G> 

600 KS0=( I -2 )*M+ 1 
Z^O.O 

IF! MO * GT # 2 ) GO TO 660 

610 KS1=KS0 
PI -DX 

Q 1 -P l * S I N2 /CO S2 
G2 = 3 * * L)Y *2 «*G 1 —P 14SIN3/COS3 
G 3-DY— CJ 1 — G2 
OO 620 JS=l*MMl 
KS1 -KS1 + 1 
Z=Z+DZ 

U ( KS 1 J^PO+Z^Pl 

620 V< KS1 >-QO + Z*<C)l*Z*( Q2+Z*Q3) ) 

GO TO 290 
660 KS2=KS0 

R=SQRT{ DX*DX+DY*DY ) 

P1=R *COS2 

P2=3.* OX-R*< 2 **COS2^COS3 ) 

P3-DX-P1 — P2 

G1=R*SIN2 

Q2=3.*DY-R *< 2.4SIN24-SIN3) 
Q3-OY-G1-Q2 
00 670 JS-1»MM1 
KS2=KS2+1 
Z=Z+DZ 

U(KS2)=P0+Z*{PH'Z*<P24Z*P3)> 
67 0 V{K$2|=00fZ<(0H2*( C2+Z*Q3J ) 

GO TO 290 

SOQ IF { MD* NE • 2 ) GC TO 890 
DO 8 10 KR=1 * N 
TS=U< KR ) 

U<KR)~V(KP> 

V(KR)=TS 
810 CONTINUE 

890 L0=LPS 
N0=N 
RETURN 
END 

SUBROUTINE SFLN2 
INTERPOLATION EY SPLINE FIT 
DIMENSION A{50)«V<£0) 
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COMMON X< 50 >fU< 5 O)iS( 5 C) .DEL < 50 ) * Q 1 . QN * N F TS 
N=NPTS 

IF ( N — 3 1 5 * 5*1 

1 IF(N-SQ) 6.6*5 

6 DEL < 2 > — X t 2 ) — X ( 1 ) 

VI 1 I (U< 2 >-U< X >>/DEL< 2 ))-Q 1 ) 

N 1 =N-I 

DO 2 1=2. N1 

DEL ( IF 1 )=X( I*t >-X( I) 

2 VtI) = i(U(I-l ) /DEL ( I ) )-U< I ) *( < 1 * 0 /DEL ( i > ) ♦( 1 * O/DEL < 1 + 1 ) ) ) 

C +(U(I + 1) / DEL C I ♦ t > ) ) * 6* 0 

V < N > = < QNt{U(Nl } - U t N ) >/CELCN> )* 6*0 
A < 1 )= 2 . 0 *DEL< 2 ) 

A(2)=l » 5*0 EL < 2 ) +2 • G* DEL < 3 ) 

V( 2 )=V( 2 )- 0 . 5 *V( 1 ) 

DO 3 1=3, Nl 
C=DEL< i)/A( 1-1 ) 

A ( I >=2 . 0 *t DEL ( UiDEL C I + 1 ) >-C*DEL ( I > 

V< I > = Vl I )-C*V( 1-1 I 

3 CONTINUE 
ODELt N ) / A < N 1 > 

A<N)= 2 . 0 *DEL( K >~C*DELt N) 

VtN)-V(N>-C*V t M ) 

StN)=V(NJ/A<N> 

DO 4 J= 1 ,N 1 
I J 

4 S(I)^(VU)-OEL(HIHS(I+l))/A{I) 

5 RETURN 
END 

FUNCTION SPLIMV) 

COMMON X( 50 >,U< 50 >*S( 50 ),DEL( 50 > .Ql.QN.NFTS 
I F C Y — X ( 1)) 1*1,2 

1 5 PLIN=U( i)-Ql*€XCll-Y> 

GO TO 5 

2 DO 4 K=2,NPTS 
IFIY-X(KI) 3,3.4 

3 CONTINUE 
K 1 — K - 1 

FFI=S(K 1 )*(X(K)-Y )**3 
FF 2 = S(K >*t Y-X< K 1 > >**3 
FF 3 = no/< 6 # 0 *DEL{K ) ) 

Fl=t FF 1 +FF 2 )*FF 3 

F 2 - IY-X<K 1 ) )*< CUf K J/DELCK 1 >-< St K>*DEL (K ) )/ 6 . 0 ) 

F3= ( X ( K ) -» Y ) ♦ { (U(K D/DELtK ))-( St K1 ) * DEL ( K) |/6.0) 

SPL I N — F 1 +F 2 +F 3 
GO TO 5 

4 CONTINUE 

SPL I N=U(NPTS ) + 0 N*( Y-X{ NPTSJ) 

5 RETURN 
END 
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APPENDIX C. CRP ERROR MESSAGES 


C . 1 Namelist Input Errors 

1. Value of variable xxxx not in correct range. Run terminates. 

Variable Correct Range 

RTYP 0,1 

ABSCA 0,1 

2. Variables DELTAU and DELTAV must both be nonzero. Run 

terminates. 

3. If linear scaling is desired, all five variables (XO, XF, 
DELTX, XFMT, and SUBDIV) must be read. Check to make sure 
XO is not equal to XF. Run terminates. 

4. Warning - graph titles have not been supplied. Blanks 
will be substituted. 


C.2 Interpolation Errors 

1. An error condition has been discovered in INTERP. The error 

return code is X, where \ 

1 = INT less than 1 or greater than M 

2 = An index in IGP(I) outside index range of 

the first through last data points of curve N. 

2. The value XI (or X2) xxx.xx does not fall on any line of 
constant X2 (or XI). 

3. The number of interpolated lines changed to NNN. The 
new array is as follows: 
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FIG . 1 


DECLINATION* EAST POSITIVE) AT SURFACE OF 
EARTH VERSUS GEODETIC LATITUDE AND LONGITUDE 
I G RF 1 9 65 MODEL , EPOCH 19 70 








Figure 4. Data Deck Setup 


Blank Card 


Bata Cards 
U=Ui 


Header Card 
U=Ui 


Data Cards \i=\J 7 


Header Card 
U=U 2 


Data Cards U=U 1 


X=f (V) for 
U-U 2 


Header Card 
U=U 1 


Blank Card 


X=f(V) for 
U=Ui 


Second Orthogonal 
Family if needed. 

RTYPE=1 only 
(Section 4.2.1) 


Data Card V=V^ 

Header Card V=Vj 

/ 


Data Cards V=V 2 


Header Card V=V 2 


X=£ QJ) for 


Data Cards V=V. 


X=f[U) fOT 
V=V, 


First Orthogonal 
Family. 


Header Card 

V=Vi 

Number of Data Cards 


ORBITAL INTEGRATION MAP EPOCH 1975 


FIG . 6 


CIRCULAR ORBIT 60 DEG INCLINATION 
SOLAR MINIMUM A E 9 AND A E 5 MODELS 



<2 






ELEC TRONS/ SQCM .DAY 



TEST OF CARPET MATRIX FOR 60 DEG. 


FIGURE 7. 



ELECTRONS/ SQCM .DAY 




TEST OF CARPET MATRIX FOR 40 K EV. 


FIGURE 8. 


O' 



ELEC TRON S/ SQCM . SEC 


TEST OF CARPET MATRIX FOR 40 KEV. 


FIGURE 9. 



ELECTRONS/ SQCM . SEC 



TEST OF CARPET MATRIX FOR 40 KE V . 


FIGURE 10. 



ELECTRONS/ SQCM . SEC 



TEST OF CARPET MATRIX FOR 100 KEV.CRVFIT 


FIGURE 


1 1 . 


-4 

<£> 



ELEC TRONS/ SQCM . SEC 





TEST OF CARPET MATRIX FOR 100 KEV. SPLINE FIGURE 12. 



ELECTRONS/ SQCM . SEC 



TEST OF CARPET F1ATRIX FOR 100 KEV. SPLINE 


FIGURE 13. 


00 



ELECTRONS/ SQCfl . SEC 



// Job Card 
// EXEC FORTRANH 
//SOURCE. SYS IN DD * 

Source Decks 

// EXEC LINKGO, REGION=170K 

//LINK. SYSLIB DD DSNAME=SYS2.SD4060,DISP=SHR 

//GO. SC4060ZZ DD DSN=job id,UNIT=7TRACK, 

// LABEL=(1,BLP) ,DISP= (NEW, KEEP), 

// DCB= (DEN=1 , TRTCH=C , RECFM=F, BLKS I ZE=1024) , 

// VOL=SER=xxxx+ 

//GO . SYSUDUMP DD SYS0UT=D 
//GO.DATA5 DD * 

Data Deck 

/* 

// 


t Tape Number. 


Figure 14. JCL Required with Source Decks 
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// Job Card 

// EXEC LINKGO, REGIONS 7 OK 

//LINK. SYSLIB DD DSNAME=SYS2.SD4060,DISP=SHR 
//LINK.SYSLIN DD * 

Object Decks 

//GO. SC4060ZZ DD DSN=job. . . id,UNIT=7TRACK, 

// LABEL=(1,BLP), DISP= (NEW, KEEP), 

// DCB= (DEN* 1,TRT.CH=C,RECFM=F,BLKSIZE= 1024) , 

// VOL=SER=xxxx^ 

//GO.SYSUDUMP DD SYSOUT=D 
//GO.DATA5 DD * 

Data Deck 

/* 

// 


t Tape Number 


Figure 15. JCL Required with Object Decks 
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LINES OF CONSTANT LLN6 


LONG - 0.0 


LEP VAR 

v 

LAT 





DEP VAR 

Y 

LAT 


3-23000 E-Qi 

10 .0 

1 .CGCCOE 

01 

4 



3. 250 00 E-01 

10.2 

1 . IOOOOE 

01 

3.27000E-0 i 

1C .4 

1 .2000011 

01 




3 • 3 0 0 0 0 E—0 1 

10.6 

1 .30 0 0 0E 

01 

3.33000E-01 

10 - J 

l .»rocoE 

01 




3.36000 E—0 1 

11.0 

i .sooooe 

01 

J.39000E-01 

11 -2 

1 . 60C00E 

01 




3 . 430 00 E—0 1 

1 1 .4 

1.70000E 

01 

3 *4 7 OOQE— 0 1 

It .6 

1 .ftOOOOE 

0 1 




3.5 00 0 0 E—0 1 

l 1.6 

1 • 90000E 

01 

3.O400C&-01 

12 .0 

2 -c ocooe 

01 




3. 5 ROD 0 E—0 1 

12.2 

2. IOOOOE 

Ol 

3.63O00E-O1 

12.4 

2.2OC00E 

01 




3.67000 E-0 l 

12.6 

2 * 30 00 OE 

01 

3 .71 000fc-0 1 

12 . a 

2 .40 CODE 

01 




3.76000E— 01 

13.0 

2.50000E 

01 

3.UOOOOE-C1 

13.2 

2.O0000E 

01 




3.6SOCOE-01 

13.4 

2.70000E 

01 

3 .69000 £-01 

13.6 

2 .tiCCOOd 

01 




3. 94000 E— Ol 

13.8 

2*90 OOCE 

01 

3.9600QE-01 

14.0 

3 .oor OOL 

0 1 




4 * 030 00 E—0 1 

1 4.2 

3. IOOOOE 

01 

4 .07000 E-01 

14.4 

J.20000c 

01 




4. 1 2000 E-01 

14.6 

3.30000E 

01 

4 . 1 6000 E—0 1 

14. b 

j.aOCOOE 

01 




4 • 2 00 00 E—0 1 

15.0 

3. SOOOOE 

Ol 

4 ,25nC0E-Ql 

lb. 2 

J.bCOOOE 

01 




4 . 290 00 E—0 1 

15.4 

3.70000E 

01 

A.33000E-01 

lb .6 

J.oOCOCE 

01 




4 « 37000 E—0 1 

15.8 

3.90000E 

01 

4.41 000fc-01 

16 .4 

4.CCOnOc 

0 1 




4 • 440 00 E—0 1 

16.2 

4. IOOOOE 

01 

4.48000E-01 

16-4 

4.20 0 00 b 

0 1 




4.520C0E-OI 

16.6 

4 .30 00 OE 

01 

4.55000t-0 1 

16.3 

4 .40 0 OPfc 

01 




4. 5 6000E-01 

17.0 

4.50000E 

01 

4.62CC0E-01 

17 .2 

4. 60 c cor. 

0 1 




4 • 650 00 E—0 1 

17.4 

4 - 700 0 OE 

01 

4 .68000E— 0 1 

IV .6 

4. '£>00 COE 

0 l 




4 . 710 00 E—0 1 

17. ti 

4.9000QE 

01 

4 . 74 000 E: -0 1 

ij .0 

b.COCCOC 

01 

LONG 

- 

40.000 

5.99000E 02 

999.0 

9 . 9900 OE 

02 

DEP VAR 

Y 

LAT 





OtP VAR 

Y 

LAT 


3 • 6200 C£— Q 1 

1W.0 

i .rcoooE 

01 




3. d5P00LL-01 

18.2 

1 .IOOOOE 

01 

3.67COOE-0 1 

13.4 

1 .20000b 

0 1 




J.7Q000E-O1 

16.6 

1 ■ 3000 OE 

01 

3.73000E-C 1 

13.0 

1 .40CC0E 

01 




3.7 60 0 0 E—0 1 

19.0 

1 . 5000 OE 

OI 

3.60000E-0 1 

lir .2 

1 . 60000 L 

01 




j.azo ooe-oi 

19.4 

1 . 7000 OE 

01 

3.87000E-01 

19.6 

1 .bOPOOE 

01 




3. y 1000 E— 01 

19.8 

1 .90OOOE 

01 

3.9S0Q0E-01 

20 .O 

2 . 00 000 E 

0 1 




3.99000E-0! 

20 .2 

2. IOOOOE 

01 

4 . 03000 E— 0 1 

20 .4 

2 . 20 C COE 

01 




4 • 070 00 E—0 1 

20 .6 

2 . 3000 OE 

01 

4.1 2 000 E— 0 1 

20 .ti 

2*4('t?00t 

01 




4. 1 6000 E— Oi 

21.0 

2 * 500 0 OE 

01 

4.20000E.-01 

21 • 2 

2.600OCE 

Cl 




4 • 240 0 0 6—01 

21 .4 

2 . 7000 OE 

01 

4 .280006-01 

21 .6 

2.r.cr-oot 

oi 




4.33000E— 01 

21 .6 

2 • 9000 OE 

01 

4.37000E-O1 

22 .O 

3 .OCCPOb 

0 1 




4.4 lOOOfc-Ol 

22.2 

3. IOOOOE 

01 

4 • 45 0 00 E—0 1 

22.4 

3. 20000b 

01 




4.490 00 E—0 1 

22.6 

3 • 30 0 0 0 E 

01 

4 .£>3000 E—0 1 

22.3 

J . 40 0 CO E 

01 




4 * b GO 00 E— 0 1 

23.0 

3.50000E 

01 

O.60000E-01 

23.2 

3 .60 CODE 

01 




4 .64500E-01 

23.4 

3.70000E 

01 

4.67000E-01 

23 .6 

J . 30 OOOL 

0 i 




4 • 70000 £—0 1 

23.6 

3.90000E 

01 

4.74000E-01 

24.0 

4 .COOOOb 

0 1 




4 . 7700 0 E -0 1 

24.2 

4. IOOOOE 

01 

4 .GOOOOE-Ol 

24 .4 

4.2C0COE 

0 1 




4 . 6 30 00 E—0 1 

24.6 

4 * 30 00 OE 

01 

4 .06OOOE-O1 

24 .6 

4 . 40C OOE 

01 




4 * 090 0 0 E—0 1 

25.0 

4. SOOOOE 

01 

4 . 92 0 0 0 E— 0 1 

25 . 2 

4 . 60 COO £ 

0 1 




4 . 950 0 0 E—0 1 

25 . 4 

4 . 700 0 OE 

Ol 

4 .90000E-01 

2b. o 

4 • JOOOOt 

01 




5. oooooe— o i 

25.6 

4.90000E 

01 

5 .03000E-01 

26 .0 

5.0000CE 

C 1 

LLNG 


30.000 

9. 99000 E 02 

999.0 

9.99000E 

02 

OLP VAR 

Y 

LAT 





DEP VAR 

Y 

LAT 


4. osoooE-a i 

26.0 

i .ooecob 

0 1 




4 .070 00 E-01 

26.2 

1. IOOOOE 

01 

4.C9000E-01 

26.4 

1 . 200 00 £ 

01 




4. 1 lOOOE-Ol 

26.6 

1 . 3000 OE 

01 

4.14000 E—0 1 

2o . 3 

1. 40000b 

01 




4. 17000 E-01 

27 .0 

l. SOOOOE 

01 

4 .2C000E-C1 

27 .2 

1 . bOOOPL 

01 




4 . 240 0 0 E—0 1 

27 .4 

1 . 7000 OE 

01 

4 • 20000E-0 1 

27 .6 

1 . 30000E 

0 l 




4.3200OE-01 

27.0 

1 • 90 0 0 0 E 

01 

4.36000E-01 

23.0 

2.000 COIZ 

0 1 




4.400 00 E—0 1 

28.2 

2. IOOOOE 

01 

4.450006-01 

20 . 4 

2 * 20 0 00 £ 

01 




4 • 500 00 E—0 1 

28.6 

2.30000E 

01 

4.S5000E-01 

2o . 3 

2.40000£ 

0 1 




4 * 60000 E—0 1 

29.0 

2. SOOOOE 

01 

4 .66000E-01 

29 ■ 2 

2.O0CCCE 

01 




4. 71000E-01 

29.4 

2.70000E 

01 


Figure 16 


Sample Line Printer Output 



Table 1. Format of Header and Data Cards 


Card 


Header 


Data 


Column 

Format Name 

Description 



1-67 

67X 

Available for Comments 


68-69 

12 N 

Number of Data Cards 

to follow 

70-74 

5X 

Blanks 



75-80 

F6.m 

Value of constant independent 



variable for following member 

1-8 

F8.5 

Independent variable 

>Pair^ 


9 

IX 

Blank 

4(n-l)+l 

10-19 

E10.3 

Dependent variable 

J 


20-27 

F8.5 

Independent variable 


4(n-l)+2 

28 

IX 

Blank 

/ Pair 

29-38 

E10.3 

Dependent variable 

) 


39-46 

F8.5 

Independent variable 

1 


47 

IX 

Blank 

} Pair 

4(n-l)+3 

48-57 

E10.3 

Dependent variable 

) 


58-65 

F8.S 

Independent variable 

| 


66 

IX 

Blank 

> Pair 

4(n-l)+4 

67-76 

E10.3 

Dependent variable 

) 


77-80 

4X 

Blanks 




Tn is the number of the data card 1 < n < N 
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Table 2. Alphabetic List of NAMELIST 
Described in Appendix A 

Parameter No. Parameter Parameter No. 


16 

ABSCA 

33 

45 

CHFMT 

38 

10 

CHOICE 

21 

43 

DCHAR 

23 

4 

DELTAU 

25 

5 

DELTAV 

39 

19 

DELTX 

11 

1 

IPR 

12 

15 

KOLOR 

8 

29 

LABE LI 

40 

34 

LABEL2 

13 

44 

NCHAR 

14 

3° 

NLAB1 

6 

35 

NLAB2 

18 

9 

NUPLOT 

46 

7 

NVPLOT 

17 

28 

NX 

31 

22 

N1 

36 

24 

N2 

41 

26 

N3 

27 

3 

PO 

32 

2 

RTYPE 

37 

20 

SUBDIV 

42 


Parameters 

Parameter 

THETA1 

THETA2 

TITLE1 

TITLE2 

TITLE3 

ULABEL 

UMAX 

UMIN 

UP LOT 

VLABEL 

VMAX 

VMIN 

VP LOT 

XF 

XFMT 

XO 

XP0S1 

XP0S2 

XTICK 

XTITLE 

YP0S1 

YP0S2 

YTICK 
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